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1.  Introduction 

The  purpose  of  this  manual  is  to  describe  the  execution  of  pre-  and  post- 
processors and  their  associated  analysis  programs  on  the  University  of  Maryland 
Univac  1103.  The  manual  contains  sections  which  describe  system  attributes 
affecting  tne  execution  of  the  programs;  and  sections  which  describe  those 
programs  availaole  in  the  library. 

The  Univac  1103  operating  under  the  EXEC  3 operating  system  supports  two 
modes  of  access  simultaneously.  These  modes  are  "Demand"  and  "Batch".  The 
demand  mode  allows  the  user  to  access  the  computer  by  means  of  a remote  terminal; 
and  with  minor  exceptions  makes  all  of  the  system  features  normally  available  to 
"Datch"  users  available  to  "demand"  users  in  a real-time  environment. 

In  the  demand  mode  tne  user  may  utilize  system  programs  for  file  management, 
for  file  editing  and  for  submission  of  programs  for  execution.  In  the  Univac  1108 
^XEC  3 environment  there  is  no  difference  between  a demand  program  and  a batch 
program;  a program  may  be  executed  in  either  mode.  However;  the  computer 
operations  staff  may  restrict  the  demand  user  to  a smaller  core  than  that  available 
in  batcn.  (The  author  has  exercised  HAST RAN  in  the  demand  mode;  and,  the  system 
responds  by  executing  the  program  immediately.) 

The  "demand"  environment  is  a new  world  for  those  of  us  who  have  graduated 
successively  from  the  earliest  programmable  computers  to  those  of  3rd  generation 
machines  with  operating  systems  which  allow  man-machine  interaction.  The  purpose 
of  developing  these  pre-processors  and  the  associated  analysis  programs  in  this 
environment  is  to  provide  a model,  however  unsophisticated,  that  will  serve  as  an 
indication  of  the  direction  that  computer  usage  is  taking. 
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2.  General  Procedures 

2.1.  Connecting  up  to  tne  machine 

The  procedure  for  connecting  to  the  computer  varies  with  the  type  of 
terminal  and  telephone  equipment  being  used.  One  first  starts  by  turning  the 
terminal  (and  sometimes  telephone  equipment)  on,  and  dialing  the  computer.  A 
busy  signal  means  that  all  incoming  lines  to  the  computer  are  in  use.  A no- 
answer may  mean  the  computer  is  shut  down,  but  during  normal  hours  usually  means 
that  the  system  has  just  "crashed"  and  is  being  re-started.  A high-pitched  "tone" 
indicates  the  computer  is  ready.  If  a teletype  with  an  integral  dial  is  being 
used,  the  connection  is  now  complete;  if  a Bell  System  data  set  is  used,  one 
snould  press  tne  "DATA"  button  and  cradle  the  handset;  users  of  acoustic  couplers 
should  place  their  telephone  nandset  in  the  coupler. 

After  establishing  the  connection  (and  before  typing  anything  else),  nit 
tne  "BREAK"  key  on  your  terminal.  If  the  system  responds  with  ^OUTPUT  INTERRUPTED*, 
you  nave  dialed  into  a telephone  line  whose  previous  user  did  not  hang  up 
properly.  Type  the  single  string  T 10'  followed  by  a carriage  return  and 

then  an  <?FIN  command  to  get  rid  of  the  previous  user's  run.  After  printing  the 
accounting  information  for  the  previous  user,  the  system  will  request  your 
account  number.  You  should  then  proceed  as  below  under  "identifying  yourself". 

If  the  system  does  not  respond  to  the  "BREAK",  proceed  to  identify  your 
terminal  as  noted  below. 
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2.1.1.  Identi fvino  vour  terminal 


r 


Once  the  connection  is  established,  you  must  identify  the  terminal  beino 
used  by  typing  in  a valid  site  I.D.  number  which  will  be  made  vailable  to 
Qualified  users.  A valid  site  I.D.  should  result  in  a response  from  the 
computer.  Occasionally,  the  system  may  not  respond  to  a valid  sit11  I.D.  If 
retvpinq  it  several  times  does  not  heln,  consult  the  "trouble  shootino"  tables 
in  Appendix  A for  corrective  action. 

Once  the  terminal  has  been  identified,  it  is  available  foe  use.  Any 
number  of  independent  ORUN's  may  be  sinned  on  and  terminated  (^FIM'ed)  while 
the  terminal  is  "active".  The  system  will  allow  a terminal  to  be  active  for  a 
brief  period  when  no  job  is  actively  using  it  (i.e.  at  initial  sion-on  or  after 
a previous  job  has  QFIN'ed);  if  no  new  DRUH's  start  within  the  allotted  time, 
the  system  will  automatically  disconnect  the  terminal. 

2.1.2.  Identifying  yourself 

Once  the  terminal  is  active,  the  system  will  ask  vou  to  identify  yourself 
by  typinq,  "ACCOUNT  NUMBER?".  In  response  to  this,  vou  mav  tvne: 

(1)  only  your  account  number 

(2)  any  or  all  fields  of  an  EXEC0.  "RUN"  card  startino  with  the 

account  number  and  continuino  with  project,  time  estimates,  °tc. 

(3)  an  entire  EXEC3  ORUN  card  image  as  shown  below: 

PRUN  <RUNID>,  <ACC0UNT  NUMBER> , <PR0JECT  ID>,  <TIME>,  <PAGES> 

where 

RUNID  is  the  six  character  identification  oiven  to  this  session. 
If  omitted  the  system  will  generate  an  I.D.  The  user  is 
cautioned  to  record  this  I.D.  for  use  in  identifying  high- 
speed printer  output. 

ACCOUNT  NUMBER  is  made  available  to  authorized  users. 

PROJECT  ID  is  a twelve  character  I.D.  which  becomes  the  file 
qualifier  if  present. 

TIME  is  the  maximum  execution  time  in  minutes. 

PAGES  is  the  maximum  number  of  pages  of  output. 


L. 


2-2 


r 


I 

'T- 

I 


In  most  cases,  alternative  (3)  is  recommended  in  order  to  provide  complete 
accounting  information.  After  entering  the  above  information,  depress  the 
carriano  return  (C-R)  to  inform  the  system  that  you  are  finished.  It  will 
respond  by  askinq  for  your  password,  which  should  then  be  tyoeu  in.  Valid 
account  numbers  and  associated  passwords  will  be  supplied  to  qualified  users. 

If  any  errors  occur  in  the  sian-on  procedure,  the  system  will  ask  for  the 
account  number  again.  Otherwise,  it  will  print  your  RU.’IID  and  today's  time  and 
date.  The  RUNID  is  internally  generated  if  options  (1)  or  (2)  are  taken  in 
enterinq  the  account  number.  If  the  entire  ORUN  card  is  supplied  then  the  RUNID 
i , taken  from  the  ORU.J  card.  The  user  is  cautioned  to  record  the  RU.'IIJ  since 
it  will  be  required  to  identify  output  that  is  directed  to  one  of  the  hiqh-speed 
ori nters . 

After  sign-on  is  complete,  the  system  will  inform  you  o*-  the  status  of 
your  workspace  (which  is  a special  file  called  the  TPF$.)  This  is  a protected 
area  of  drum  which  the  system  establishes  for  you  as  a work  area,  .’formally 
the  workspace  is  discarded  when  you  terminate  a terminal  session  and  a new  one 
is  established  at  the  next  terminal  session.  If  the  system  crashes  durino  a 
terminal  session,  however,  the  workspace  will  be  retained  and  recovered  the 
first  time  you  sign  on  after  the  crash. 

Since  system  crashes  rarely  involve  "down"  times  of  more  than  a few 
minutes,  workspaces  which  are  saved  because  of  a system  crash  are  short-lived. 
They  will  usually  be  immediately  available  if  one  signs  back  on  to  the  machine 
on  the  day  of  the  crash.  If  not  recovered  on  that  dav,  they  will  be  rolled  out 
onto  tape,  and  a short  wait  will  be  required  to  recover  then.  Workspaces  not 
recovered  within  tv/o  days  after  being  rolled  out  to  tape  are  nurqed  from  the 
system. 
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2.2.  Tyninq  Terminal  Input 
2.2.1.  Editing  Input  Lines 


Input  to  the  system  from  a demand  terminal  is  nrocessed  a line  at  a tine  as 
it  is  tyoed;  each  line  must  be  followed  by  a C-H  before  it  will  be  accepted.  At 
any  time  before  the  C-R  key  is  typed,  the  line  may  bo  edited.  Tyoinn  a question 
mark  will  cause  the  entire  line  to  be  deleted;  tvoinq  an  underline  (back  arrow 
on  some  teletype  terminals)  will  cause  deletion  of  the  last  non-deleted 
character  on  a line  (thus  the  last  four  characters  are  deleted  by  tyninq  four 
underlines).  Although  the  system  will  accept  either  unner  or  lower  case  i inut, 
lower  case  letters  will  be  translated  into  unner  case. 

2.2.2.  ASCII  Control  Codes 

When  in  normal  (not  full  ASCII)  input  mode,  many  of  the  ASCII  communications 
control  codes  have  special  meaning.  On  teletype  terminals,  these  codes  are 
usually  generated  by  holding  down  the  CTRL  key  and  striking  an  alphabetic  key, 
although  some  (like  ESC)  have  a key  of  their  own.  Control  codes  not  noted  in 
the  table  below  are  translated  by  the  system  into  question  marks.  The  codes 
which  are  used  are: 


CODE 

KEY 

USE 

XON 

CTRL-Q 

start  form  I naper  tane  inout 
(see  EXEC3  PRM) 

XOFF 

CTRL-S 

terminate  form  I mode  (must  o 
sent  twice) 

TAPE 

CTRL-R 

start  form  II  oaner  tape  inout 
(see  PRM) 

T-A-P-E 

CTRL-T 

terminate  form  II 

SHIFT-OUT 

CTRL-M 

start  full  ASCII  input  (see  UOM 
documentation) 

SHIFT - IN 

CTRL-0 

terminate  full  ASCII  input 

I 
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m or  ENQ 

CTRL-E 

causes  system  to  echo  its  sign-on 
line,  useful  in  determining  whether 
system  is  "up" 

EOT 

CTRL-D 

disconnects  terminal  from  computer 

RUBOUT 

RUBOUT 

ignored 

BREAK 

BREAK 

interrupts  currently  executinq  nroqram, 
see  below 

ESC 

CSC 

can  be  used  before  typino  question 
mark  to  allow  input  of  question  mark 
rather  than  deletion  of  line 

2.2.3.  Interrupting  the  Computer 

The  BREAK  key  on  teletype  terminals  is  used  to  interrupt  the  computer . When 
BREAK  is  hit,  the  computer  will  respond  with  ‘OUTPUT  INTERRUPTED*,  temporarily 
suspend  what  it  is  doing,  and  wait  for  inout  describing  further  action  to  be 
taken.  The  currently  defined  actions  are: 

@@X  JI0  Abort  whatever  is  goinq  on  and  return  control 

immediately  to  the  system  executive.  If  any 
printout  was  produced  prior  to  the  break,  it 
will  be  printed  before  the  system  waits  for 
input. 


? Ignore  the  Break  and  continue.  Especially 

useful  when  using  CRT  display  terminals  to 
"freeze"  the  contents  of  the  screen  for 
examination  before  proceeding. 

C-R  Interrupt  the  current  action  of  the  executino 

program  and  wait  for  innut.  This  action  is 
specific  to  the  9ED,  9APL,  and  OBASIC 
processors.  In  other  circumstances,  BREAK-C-R 
is  equivalent  to  BREAK-X. 
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2.3.  System  Functions 


In  order  to  utilise  the  terminal  device  to  access  the  computer  and  execute 
tne  program  the  user  must  have  the  capability  to: 

(1)  Create  files 

(2)  Move  files 

(3)  Purge  files 

(4)  Edit  files 

(5)  Call  for  proqram  execution 

The  treatment  qiven  to  these  commands  will  be  sufficient  to  allow  the  use 
of  the  system,  but  will  be  brief  and  incomplete.  The  user  should  consult 
references  (3)  and  (4)  for  more  detail. 

2.3.1.  Files  and  elements 

Users  familiar  with  other  timesharino  systems  should  note  that  the  EXEC3 
concent  of  a "file"  is  more  aoneralized  than  that  of  many  other  systems.  A 
file  is  a collection  of  entities  called  elements.  Since  files  contain  elements 
the  file-naming  procedure  must  be  capable  of  identifying  a specific  element  of 
a file.  The  file  name  convention  is  as  follows: 

QUAL I F I ER  *F I L ENAME . EL  EMENTNAME 

The  various  parts  of  the  name  may  be  omitted  by  the  user  with  default 
values  as  follows. 

QUALIFIER  Taken  to  be: 

(1)  The  project  1.0.  if  giver  on  the  3RUN  card 

(2)  The  first  six  digits  of  the  account  number 

FILENAME  Taken  to  be  a special  file  called  the  TPF$  if  omitted. 

The  TPF$  is  the  user's  workspace. 
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Suppose,  for  example,  that  a user  had  logged  on  using  the  project  I.L. 
SYMPOSIUM",  and  wished  to  identify  the  element  ELMT  in  the  file  SYMPOSIUM* 


TPF$.  ; then  the  specification  of  ELMT  in  a suitable  system  command  would  oe 
sufficient,  since  the  system  would  supply  the  correct  default  names  for  the 
qualifier  and  file  name. 

The  user  should  note  that  the  file  is  the  basic  cataloged  entity  on  tne 
EXEC8  system.  A file  is  almost  always  indicated  by  a terminating  period  (.). 

The  period  may  be  omitted  only  when  the  system  would  understand  that  a file  is 
designated  if  omitted. 

2.3.2.  Creating  and  Cataloging  Files 

The  following  sequence  of  commands  will  catalog  (save)  a file: 

@CAT ,J  FILE , F 
(SFREE  FILE. 

0ASG.A  FILE. 

At  this  point  FILE  exists  but  is  empty.  Data  or  elements  may  now  be 
entered  into  the  file  by  using  the  text  editor  a described  in  section  2.3.6. 

2.3.3.  Moving  Files 

The  user  may  operate  directly  on  a file  in  permanent  storage,  or 
alternatively  on  a copy  of  the  file  in  the  TPF$.  The  @C0PY  command  is  used 
to  copy  the  file's  contents  from  one  storage  location  to  another.  The  form  of 
command  is: 

(SCOPY  <SPEC1>,  <SPEC2> 
where 

SPEC1  is  the  file  name  of  the  file  contents  to  be  moved,  if  blank 
the  TPF$  is  assumed. 

SPEC2  is  the  file  name  for  the  destination  of  the  file  contents, 
if  omitted,  the  TPF$.  is  assumed. 
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2.3.4.  Assigning  Files 


r 


Files  may  be  created  for  the  duration  of  the  run  by  the  command 
3ASG.T  FILENAME, F 

A file  tnat  already  exists  may  be  assigned  to  the  present  run  by  the 

command 


PASG.A  FILENAME. 


2.3.S.  Purging  Files 

The  user  may  purge  a permanent  cataloged  file  by  the  command 
0OELETE  FILENAME. 


If  the  file  is  temporarily  assigned  to  the  run,  it  may  be  purged  by  the 

conmand 


(JFREE  FILENAME. 

The  PFREE  command  frees  a file  from  the  current  run.  If  it  is  desired  to 
erase  the  current  contents  of  a file  the  following  command  may  be  used 
PERS  FILENAME. 

2.3.6.  Text  Editor 

The  University  of  Maryland  text  editor  Is  described  by  reference  2.  Ne  will 
give  a orief  description  of  the  features  which  are  required  to  create  and  edit 
data  files. 

The  editor  is  called  by  the  command  PED  <Spec>,  where  <$pec>  refers  to  an 
element  or  a file.  The  editor  operates  in  two  modes;  EDIT  and  INPUT.  If  the 
file  is  empty,  the  editor  assumes  the  INPUT  mode;  otherwise,  the  EDIT  mode  is 
assumed.  Mode  change  is  affected  by  a C-R. 

2. 3.6.1.  Input  mode 

In  this  mode  any  character  date  entered  on  tne  keyboard  is  entered  in  the 
file.  Each  physical  line  of  input  is  terminated  by  a C-R.  The  system  editing 


r 


T 


features  described  in  section  2.2.1.  are  available  in  this  mode.  The  input  mode 
is  terminated  two  successive  C-R's;  if  an  3-sign  is  entered  in  either  mode  the 
editing  session  is  terminated  and  control  is  returned  to  the  system  executive. 

2. 3. 6. 2 Edit  ilode 

In  this  mode  the  editor  prompts  the  user  with  an  asterisk  (*)  to  enter  an 
edit  command.  The  U of  M text  editor  has  a great  deal  of  sophistication; 
however  a few  commands  will  suffice  for  our  purposes.  The  form  of  the  command 
is  a command  mnemonic  followed  by  a blank  with  control  information. 


Command 
G LN 

L String 
N L 
U L 


Description 

Go  to  and  print  line  number  LN 

Locate  and  print  the  first  line  in  which  "string"  occurs 
Move  down  L lines  and  print 
Move  up  L lines  and  print 


C #String#Stringl#  At  the  current  line,  change  String  to  Stringl. 

The  delimiters  (#)  may  be  an^  character  not 
included  in  the  string. 


P L 
0 L 
0 L 
T . 
E 


Print  L lines  of  text  starting  with  the  current  line. 
Print  L lines  of  text  starting  with  the  next  line. 
Delete  L lines  of  text,  starting  with  the  current  line. 
Position  the  pointer  to  the  top  of  the  file 
End  the  editing  session 


2.3.7.  Redirecting  Output 

The  program  output  in  either  mode  may  be  redirected  (from  the  terminal  or 
tne  line  printer  for  demand  and  batch,  respectively)  by  the  use  of  the  3BRKPT 
command.  The  form  of  the  command  is, 

3BRKPT  PRINTS/FILENAME 
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where  FILENAME  is  NOT  terminated  by  a peri od. 

The  file  must  be  cataloged. 

If  we  desired  to  redirect  tie  output  to  a file^called  OUT,  for  example, 
tne  following  sequence  of  commands  could  be  used: 

^CAT.J  0UT,F  (catalog  the  file) 

i$FREE  OUT.  (free  tne  file) 

@ASG,A  OUT.  (assign  with  A-option) 

3BRKPT  BRINTS/OUT  (direct  the  print  file  to  OUT) 

On  the  teletype  al_[  output  is  now  redirected  to  OUT.  To  redirect  tie 
output  back  to  the  teletype  type: 

$BRi<?T  PRINTS 


2.3.3.  Program  Execution 

Tne  EXEC  3 features  which  are  available  to  tne  demand  user  are  discussed  in 
references  (5)  and  (4).  We  will  present  only  those  commands  which  are  re qo  red 
for  the  execution  of  programs  in  either  batch  or  demand  mode.  The  basic 
difference  in  the  two  modes  is  that  demand  programs  are  executed  immediately 
in  real  time  while  batch  requests  are  submitted  to  the  batch-queue.  The  advantage 
of  the  latter  is  that  the  terminal  is  free  to  perform  other  tasks  during  batch 
execution  of  the  program. 

2. 3. 8.1.  demand  Execution 

The  symposium  input  data  processors  (IDP)  and  their  associated  analysis 
programs  are  saved  on  the  files 

3YMP0S I UM*P  REPROCESSOR . 


and 

SYMP0SIUM*PR0GRAM. 
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The  user  may  interrogate  either  file  to  determine  the  availaole  executable 


modules  oy  using  tne  system  contnand. 
tfPRT.T  < F I LENAME  > 

After  determining  the  contents  of  the  files,  any  of  the  absolute  modules 
may  be  executed  directly  in  the  demand  mode  by  the  system  command 
cJXQT  <SPEC> 

where  <SPEC>  is  the  complete  element  name  for  the  executable  module.  For  example 
the  element  SAP  exists  in  the  preprocessor  file,  and  may  be  executed  by  the 
command. 

OXQT  SYMPOSIUM*PREPROCESSOR. SAP 
2. 3. 8. 2.  Batch  Execution 

There  are  at  least  two  ways  to  submit  a job  to  the  batch  queue  for  execution; 
these  are  the  0BATCH  and  3START  system  commands. 

The  $BATCH  comnand  has  the  following  form. 

3BATCH  <SPEC> , <DATA> , RMENGR 

where  *SPEC>  is  the  absolute  module  to  be  executed,  and  <DATA>  is  the  name  of  the 
cataloged  data  file  containing  the  input  data.  This  command  is  relatively  simple; 
but,  it  will  cause  the  output  of  the  analysis  program  to  be  printed  on  one  of  the 
system's  higher  speed  printers.  Since  there  is  no  means  of  redirecting  the 
program  output  to  a data  file  this  command  should  be  used  only  by  those  users 
who  have  ready  access  to  University  of  Maryland  physical  facilities. 

The  START-  command  on  the  other  hand  allows  the  user  more  flexibility;  but, 
at  the  price  of  having  to  construct  a file  containing  the  run  stream  images 
required  to  perform  tne  required  task.  The  form  of  the  command  is: 

(?START  <spec> 
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where  <spec>  is  the  name  of  an  element  or  a file  containing  run  stream  images. 
The  START-  command  simply  inserts  a file  containing  run  stream  images  in  the 
batch  queue;  the  user  must  define  the  run  stream  images  required  to  control  the 
computer  data  management  and  execution. 

For  example,  assume  that  we  had  generated  a data  file  "INPUT *'  using  the 
IDP  for  the  analysis  program  SAPIII.  We  would  like  to  execute  the  SAPIII 
program  and  redirect  the  output  to  the  file  OUT  which  you  could  inspect  at  some 
later  time  using  the  text  editor. 

We  would  begin  by  first  cataloging  the  file  which  will  contain  the  run- 
stream  images;  let  us  call  this  file  RUNIT.,  then  we  input  the  following  system 
commands: 

3CAT.U  RUN IT 
3FREE,  RUN IT. 

@ASG,A  RUN IT. 

Next,  we  will  use  the  editor  to  enter  the  run  stream  images  into  RUNIT: 
@ED  RUNIT. 

#RUN  <RUNID> , <ACCT  ,‘J0>,  <PR0JECT  I.D.>,  <TIME>,  <Pages> 

#A$G,U  0UT  (catalog  the  output  file) 

#FREE  0UT. 

#ASG,A  0UT. 

#BRKPT  PRIMTS/0UT  (redirect  the  output) 

#XQT  SYMPOS I UM*PR0GRA i . SAP III  (execute  SAPIII) 

#ADD  <qualifier>*INPUT.  (add  the  data  file) 

#FIN  (required  to  terminate  this  job) 

(C-R)  (carriage-return  to  enter  edit  mode) 


We  would  now  be  in  the  "EDIT"-mode.  Note  that  each  of  the  control  state- 
ments has  been  preceeded  by  a #-Sign  rather  than  an  @-sign.  This  was  done 
because  the  3-sign  would  have  caused  the  control  to  return  to  the  system.  We 
will  therefore  change  the  # to  an  Q-sign  and  list  the  file  at  the  same  time  by 
the  editing  commands; 

T - go  to  the  top  of  the  file 

c»*  /#/3/  _ change  the  # to  @ in  every  line  in  the  file 

To  exit  the  editor  we  now  enter 

E - end  the  edit  session 

Tne  required  suDmission  can  now  be  accomplished  by  the  system  command 
3START  RIJNIT. 

2.3.9.  Printing  File  Contents 

A hard  copy  of  the  file  contents  can  be  obtained  by  using  a system  command 
to  request  that  it  be  printed  on  a high-speed  printer.  The  command  is: 

3SYM  <SPEC>,,  RMEUGR 

The  file  "OUT"  created  by  the  analysis  program  SAP III  could  be  printed 
by  the  command 

@SYM  OUT.,,  RMENGR 
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3.  Program  Function 
3.1.  Input  Data  Processors 
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The  input  data  preprocessors  (IDP)  are  computer  programs  that  are  normally 
executed  in  the  demand  mode  and  which  allow  the  user  to  create  an  input  data  file 
for  the  appropriate  analysis  program.  The  preprocessors  operate  in  either  of 
two  modes;  prompting  or  data  file. 

In  the  prompting  mode  the  IDP  carries  on  a dialogue  with  the  user;  it 
asks  for  input  data  sets,  and  defines  the  input  data  items.  The  IDP  then  writes 
a formated  output  file  which  is  the  input  file  for  the  appropriate  analysis  program. 

In  tne  data  file  mode  the  user  preprepares  the  responses  to  the  IDP 
prompting  requests  in  a data  file  or  data  element  using  the  text  editor  (see 
Section  2.3.6).  The  user  then  executes  the  IDP  and  in  response  to  the  request, 

"Do  you  want  to  enter  data  from  a data  fi le",  responds  "yes".  The  IDP  will  then 
request  tne  name  of  the  data  file  or  element  that  contains  the  input  data  file. 

The  IDP  input  will  then  be  read  from  this  file  without  the  prompting  that  is 
associated  with  the  interactive  mode.  The  system  may  be  directed  to  read  the 
responses  to  the  IDP  prompting  from  a precreated  data  file  by  using  the  command 
CjADD  <spec> 

where  <spec>  is  the  name  of  the  data  file  or  element. 

The  IDP's  presented  oy  this  symposium  generally  lack  the  capability  of 
recovering  from  a data  transmission  error.  Such  an  error,  which  may  be  due  to 
a data  generation  error  or  problem  in  long  distance  communications  lines,  will 
cause  a normal  exit  from  the  IDP.  The  only  option  available  to  the  user  is  that 
of  re-executing  the  program.  Prior  to  re-executing,  the  data  file  that  was 
created  in  the  previous  execution  should  be  deleted  by  means  of  the  system  command: 
^DELETE  FILENAME. 


L 
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The  use  of  the  prompting  mode  seems  most  desirable  as:  a teaching 

device  to  help  the  user  determine  the  sequence  and  type  of  data  required;  as  a 
means  of  inputting  unformatted  input;  and,  as  a means  of  checking  data  prior 
to  execution  of  the  analysis  program.  The  more  normal  usage  would  seem  to  De 
that  of  the  data  file  mode.  If  the  attempted  execution  of  the  IDP  in  the  data 
file  mode  leads  to  an  undecipherable  error  message  it  is  suggested  that  the 
prompting  mode  be  utilized  with  the  data  file  3ADD'ed  to  the  run  stream.  This 
will  cause  the  IDP  to  print  all  prompting  but  to  read  the  input  from  the  data 
file.  This  may  assist  the  user  in  determining  the  incorrect  entries  in  the 
data  file. 

3.2.  Analysis  Programs 

The  analysis  programs  are  executable  in  either  the  batch  or  demand 
environments;  subject  to  site  core  size  limitations.  The  input  data  for  the 

programs  may  be  prepared  by  means  of  the  appropriate  IDP;  or  in  the  case  of 

a program  wnich  requires  a formatted  input  file  and  for  which  no  IDP  has  been 
created,  by  means  of  the  text  editor. 

If  the  program  is  executed  in  the  demand  mode  the  data  file  is  added  to 

the  input  data  stream  by  means  of  an  @ADD  <SPEC>  - command.  This  command 

causes  the  EXEC8  system  to  use  the  <spec>-  file  for  input  rather  than  the  normal 
system  input  device. 

The  analysis  program  output  will  normally  be  returned  to  the  terminal 
device.  In  the  case  of  one  of  the  small  special  purpose  analysis  programs  this 
may  be  acceptable.  However,  the  output  associated  with  most  of  the  larger  programs 
is  quite  extensive  so  that  some  means  of  redirecting  the  program  output  is  desirable. 

This  feature  called  "Breakpointing"  is  discussed  in  section  2.3.7.  You  should 
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carefully  note  that  once  the  output  has  been  redirected, all  computer  responses 
will  be  written  on  that  file  until  the  system  has  been  informed  to  do  otherwise. 
The  command  iaBRKPT  PRINT$  must  be  given  to  direct  the  output  back  to  the  terminal. 

The  analysis  program  may  be  sufficiently  large  that  it  cannot  be  executed 
in  the  demand  mode;  or,  requires  more  execution  time  than  is  reasonable  for 
demand  submission.  In  this  case  the  job  should  be  submitted  to  the  batch  queue 
using  tne  system  commands  discussed  in  section  2. 3. 8. 2 
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4.0  SAP  - Structural  Analysis  Program 

4.1  Introduction 

SAP  III,  A Structural  Analysis  Program  for  Static  and  Dynamic  Response  of 
Linear  Systems  is  described  in  reference  (4).  Tne  analysis  program  is  executed 
on  the  Univac  1108  under  control  of  the  EXEC3  operating  system  using  single 
precision  arithmetic. 

Tne  program  has  the  following  analysis  capability: 

1.  Static  Analysis 

2.  Natural  Modes  and  Frequencies 

3.  Modal  Response  Spectrum 

4.  Modal  Dynamic  Response 

Tne  structural  systems  to  be  analyzed  may  be  composed  of  comoinations  of 
a numoer  of  different  structural  elements.  The  program  presently  contains  the 
following  element  types: 

1.  three-dimensional  truss  element 

2.  three-dimensional  beam  element 

3.  plane  stress  and  plane  strain  element 

4.  two-dimensional  axisymmetric  solid 

5.  three-dimensional  solid 

6.  thick  shell  element 

7.  thin  plate  or  thin  shell  element 

8.  boundary  element 

9.  pipe  element  (tangent  and  bend) 

The  program  capabilities  and  theoretical  oasis  are  described  in  reference  (4). 
The  purpose  of  this  section  of  the  report  is  to  describe  (1)  the  IDP  which 
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generates  the  formated  input  file  for  SAPIII;  and  (2)  execution  of  tne  SAP III 
UP  and  analysis  programs  on  the  University  of  Maryland  Univac  1103. 

4.2  SAPIII  Input  Data  Preprocessor 

Tne  SAPIII  IDP  operates  in  the  Demand  mode  and  allows  tne  user  to  supply 
unformated  input  data  in  a program-controlled  logical  sequence.  The  IDP  then 
performs  some  elementary  checks  and  writes  a file  which  contains  the  formated 
SAPIII  input.  The  IUP  executes  in  either  of  two  modes,  interactive  or  data  file. 

In  the  interactive  mode  the  user  supplies  tne  required  data  in  response  to 
program-generated  prompting.  In  the  data  file  mode  the  user  prepares  a file 
containing  tne  UP  input  in  exactly  tne  same  sequence  as  that  required  by  tne 
interactive  mode.  The  IDP  then  by-passes  all  interactive  responses  and  reads 
tne  data  from  the  pre-prepared  file.  In  either  input  rode  the  preprocessor  will 
print  a sumnary  of  tne  SAPIII  data.  The  sumnary  data  is  obtained  by  reading  tne 
FORMAT’ed  SAPIII  input  file  created  by  the  IDP;  and,  thus  provides  a means  of 
checking  the  correctness  of  tnis  file. 

Tne  UP  logically  omits  requesting  redundant  input  data;  and,  logically 
will  not  print  redundent  information  in  the  output  summary.  Thus,  if  the  user 
has  specified  the  static  option,  the  UP  will  not  request  parameters  which 
pertain  only  to  the  dynamic  analysis.  The  IDP  will  insert  all  required  data 
including  blank  card  images  in  the  SAPIII  input  file. 

The  input  data  items  required  by  the  SAP  preprocessor  are  shown  by  Table  4.1. 
This  data  may  be  entered  under  IDP  in  the  interactive  mode;  or  a data  file  may  be 
created  and  referenced  uy  appropriate  response  to  IDP  interrogation  as  descr  bed 
by  section  3. 1 . 
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4.2.1.  Interactive  mode 

In  the  interactive  mode  the  IDP  describes  the  input  data  requirements  so 
that  the  experienced  user  may  describe  a physical  problem  in  terms  of  SAPIII 
parameters  with  little  or  no  reference  to  the  SAPIII  user's  manual.  It  is 
highly  recommended  that  the  inexperienced  analyst  use  the  interactive  mode  in 
conjunction  with  the  SAPIII  manual.  The  user  will  note  that  the  IDP  generally 
has  the  same  variable  names  as  those  given  in  the  SAPIII  manual;  and,  that  the 
logical  sequence  of  input  is  the  same. 

In  either  mode  the  IDP  requests  the  name  that  the  user  wishes  to  assign  to 
the  IDP  output  file.  Upon  termination  of  the  IDP  the  named  file  will  contain 
the  formatted  input  data  file  for  the  SAPIII  analysis  program.  This  file 
is  cataloged  under  the  users'  account  number  as  a private  file.  The  non-1108 
user  should  note  that  the  entire  file  name  is: 

QUALI FIER*MAME 

Where  the  qualifier  is  that  which  was  generated  at  sign-on  time. 

The  user  may  enter  the  file  name  in  full  including  any  qualifier  and 
read  and  write  keys  as  follows: 

QUAL I F IER*NAME/KEY 1/KEY2 
Where  KEY1  is  a write  key 
KEY2  is  a read  key 

The  user  is  cautioned  to  retain  these  keys  in  a secure  place,  if  provided; 
especially  if  a CRT  terminal  is  used. 

4.2.2.  Data  File  (non-prompting  mode) 

In  this  mode  the  analyst  prepares  a data  file  containing  the  data  described 
by  Table  4.1.  The  elements  of  the  data  file  are  read  by  the  IDP  in  the  same 
logical  sequence  used  by  the  interactive  mode.  In  the  data  file  mode  the 
prompting  is  bypassed. 
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4.2.3.  Description  of  IDP  Input  Requirements 

Tne  data  elements  defined  on  Table  4.1  are  generally  given  the  same  names 
as  tnose  found  in  the  SAP  users  manual;  and,  the  user  is  directed  to  eitner  the 
users  manual,  reference  (4),  or  to  the  interactive  execution  of  the  program  for 
their  description. 

There  are  some  differences  between  the  SAPIII  batch  program  input  require- 
ments, and  tne  data  set  requested  by  the  IDP.  These  are: 

1.  Data  which  applies  to  the  dynamic  analysis  are  not  requested 
if  the  number  of  loads  (LL)  is  greater  than  zero. 

2.  Data  which  describes  element  loads  are  not  requested  if  the 
dynamic  option  (LL.EQ.O)  is  chosen. 

3.  Element  loads  can  be  set  to  zero  by  answering  "yes"  to  an 
appropriate  program  query. 

4.  .lode  data  input  terminates  when  the  node  number  equals  the 
total  number  of  nodes,  NUMNP. 

5.  Element  data  terminates  when  the  element  number  equals  the 
total  number  of  elements  specified  in  a group. 

6.  The  user  may  specify  an  isotropic  material  for  the  shell 
element;  the  preprocessor  then  calculates  the  material 
stiffness  matrix. 

7.  All  data  must  be  entered  correct  as  to  type;  i.e.  real  numbers 
must  have  decimal  point.  The  normal  FORTRAN  convention  for 
variable  types  is  followed,  except  as  noted  in  the  variable 
description  which  is  given  by  the  interactive  mode. 

8.  Numerical  codes  are  used  in  place  of  holorith  data  for 

(a)  Node  point  coordinate  system  specification. 

(b)  Connect  and  Load  data  for  pipe  element. 
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4.2.4.  Example  - Uniformly  Loaded  Cantilever  Beam 


r 


Consider  the  problem  of  a beam  subject  to  a uniform  load  as  shown  below 


■\ — ! — 1 — 1 — 1 I i l — 1 I l — I 1 

/O* 


£=  n 

/r  — / • /. 

Z = ///■ 


We  chose  to  descretize  tne  beam  as  shown. 
j 4 


Because  of  the  structural  topology  we  can  use  both  node  point  and  element 
generation  features  of  the  SAP-program.  Node  point  data  for  nodes  2 through  6 
will  be  automatically  generated  by  the  analysis  program  if  on  the  data  associated 
with  node  2,KN=0  and  for  node  6,  KN=1 . 

Similarly  the  element  data  for  elements  1 through  5 may  be  generated  by 
setting  KNal  on  the  input  data  for  element  1. 

4.2.4. 1 Interactive  Mode 

Before  calling  for  execution  of  the  preprocessor  we  note  that: 

(1)  NUMNP  = 7 

NETYP  a 1 

LL  = 1 

MODEX  = 0 

KEQB  = 0 

(2)  Rectangular  Coords  will  be  used 

(3)  That  node  1 and  7 are  fully  fixed; 

that  nodes  2 through  6 must  have  degrees  of  freedom 
numbers  1,  3,  4 and  5 constrained  in  order  to  make  the 
stiffness  matrix  nonsingular 

that  node  7 is  required  to  define  the  local  X-Y  plane 


L. 
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(4)  Tne  distributed  load  must  be  modeled  using  either  fixed 
end  forces  or  concentrated  nodal  forces  - we  choose  to 
use  the  fixed  end  forces 

The  interactive  preprocessor  execution  is  as  follows,  where  the  user 
responses  to  program  or  system  requests  are  underlined. 


0XQT  SYMPOS IUM*PREPROCESSOH . SAP 


*****  TS0006***** 

INPUT  DATA  PREPROCESSOR  FDR  S.A.P. 

UNIV/AC  1100-SERIES  VERSION.  INPUT  DATA  IS  TD  BE  ENTERED  FRIT'i  EITHER 
THE  INTERACTIVE  OR  THE  DATA  FILE  MODE  USING  FREE  FORM  INPUT  AS 
DESCRIBED  IN  THE  UNIVAC  FORTRAN  V MANUAL.  IN  THIS  MTDE  THE 
INDIVIDUAL  ELEMENTS  OF  THE  DATA  SET  ARE  SEPARATED  BY  COMMAS  <»>• 
NOTE  THAT  THE  DECIMAL  POINT  <•»  IS  REQUIRED  FOR  REAL  NUMBERS.  THh 
NORMAL  FORTRAN  CONVENTION  IS  FOLLOWED  FOR  VARIABLE  NAMES  UNLESS 
OTHERWISE  INDICATED  BY  APPENDING  AN  ANNOTATION  TO  THE  VARIABLE 
DESCRIPTION 


ENTER  NAME  OF  OUTPUT  DATA  FILE 
SAP  IN 


DO  YOU  WANT  TO  ENTER  DATA  FROM  A DATA  FILE  (YES  OR  NO) 
NO 


ENTER  TITLE  (NOT  TO  EXCEED  72  CHARACTERS) 
UNIFORMLY  LOADED  CANTILEVERED  BEAM 


-MASTER  CONTROL  INFORMAT  I ON- 

ENTE*™E  SET  (NUMNP.NELTYP.LL.MODEX.KEQB)  WHERE  I 

NUMNP  -NUMBER  OF  NODES  M 

NELTYP- NUMBER  OF  ELEMENT  TYPES 
LL  NUMBER  OF  STRUCTURAL  LOADS 

• GE.l  STATIC  LOAD  CASE 

* EQ.O  DYNAMIC  LOAD  CASE 
MODEX  PROGRAM  EXECUTION  MODE 

•EQ.O  PROBLEM  SOLUTION 
•EQ.l  DATA  CHECK  ONLY 
KEQB  SET  EQUAL  TO  ZERO  (SEE  NOTE  8) 

7.1. I. O.Q 
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-NODAL  POINT  DATA 


ENTER  SUFFICIENT  DATA  SETS  OF  THE  FOLLOWING  FORM  TO  DEFINE  7 NODES 
C IC»N«  < IXCN* I ># 1=1 *6)»X*Y»Z> KN»  T ) WHERE* 

IC  IS  A COORD  CODE 

• EQ  • 1 RECTANGULAR 

• EQ.2  CYLINDRICAL 
N =NODE  NUMBER 

IX-F  IX  ITY  CODE  FOR  D.  0»  F (I> 

• EQ • 0 NO  CONSTRAINT  ON  D.p.F  Cl) 

• EQ • 1 D.O.F  (I)  CONSTRAINED 

CX»  Y*  Z ) COORDS  OF  NODE  N 
KN-NODE  DATA  GENERATION  CODE 
T -TEMPERATURE 


-ELEMENT  DATA- 

ELEMENT  DATA  MUST  BE  SPECIFIED  FOR  1 ELEMENT  GROUPS 


ENTER  DATA  SET  < NTYPE. NELMT. NMAT > FOR  GROUP  1 WHERE: 
NTYPE-  ELEMENT  TYPE 
NELMT* NO  OF  ELMTS  IN  GROUP 
NMAT  -NO  OF  MATERIAL  PROPERTIES 

1.5.1 


-THREE  DIMENSIONAL  BEAM  ELEMENT- 

ENTER  THE  DATA  SET  CNPROP»NFIX>  WHERE* 
NPR OP -NUMBER  OF  PROPERTY  SETS 
NF IX  -NUMBER  OF  FIXED-END  FORCES 

LlO 


ENTER  1 MATERIAL  CARDS  OF  THE  FORM  ( M I D» E, NU» RHO# WT > WHERE* 
MID  -MATERIAL  I.D. 

E -MODULUS  OF  ELASTICITY 
NU  -POISSONS  RATIO 
RHO  -MASS  DENSITY 
WT  -WEIGHT  DENSITY 
. 1 . E7 . >3.Q«.Q. 
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ENTER  1 PROPERTY  CARDS  OF  THE  FORM  ( P I D#  AREA*  SHR2# SHR3# IT OR# I?»  13) 
WHERE! 

PID  -PROPERTY  I • D.  (INTEGER) 

AREA=AX IAL  AREA 

SHR2-SHEAR  AREA  IN  LOCAL  2-DIRLCTION 
SHR3-SHEAR  AREA  IN  LOCAL  3-DIRECTION 
I TOR  = T OHS  I ONAL  CONSTANT  (REAL) 

12  -FLEXURAL  INERTIA  ABOUT  2-DIRECTION  (REAL) 

13  -FLEXURAL  INERTIA  ABOUT  3-DIRECTIfN  (REAL) 

1 * 1 .*Q.*0.*  I .*  1 .*  1 . 


DO  YOU  WANT  TO  SET  ALL  ELEMENT  LOAD  FACTORS  TO  ZERO  (YES  nR  NO) 
YES 


BEAM  CONNECTIVITY 

ENTER  SUFFICIENT  DATA  SETS  ( NELM* I » J. K. M I D» P I D# I DA# I DB# I DC # I DD» I REL , . IRE 
L* KG  ) 

TO  DEFINE  5 ELEMENTS  WHERE! 

NELM-ELEMENT  NO* 

( I # J*  K >- NODE  NUMBERS 
MID  -MATERIAL  I*D* 

PID  -PROPERTY  I.D.  (INTEGER) 

( IDA#  I Dll#  I DC#  I DD )-  FIXED-  END  FORCE  ID  S FOR  LOAD  CASES  (A.B.C.D) 

( I REL # JR EL )- END  RELEASE  CODE 
KG  -ELEMENT  GENERATION  PARAMETER 
NOTE!!*  DATA  SET  FOR  LAST  ELEMENT  MUST  BE  DEFINED 
I »I*2*7»1*1»0#Q#0*0»0*0*1 
5*5»6#7»1#1*0*0*0*0*0*0*Q 


-CONCENTRATED  LOAD /MASS  DATA 

ENTER  THE  DATA  SET  ( N* L* FX# FY» FZ»MX# MY»MZ > WHERE! 

N -NODE  NO. 

L -STRUCTURAL  LCttD  CASE 

( FX»  FY»  FZ ) - FORCES  OR  TRANSLATIONAL  MASSES 
( MX  * MY  # MZ ) - MCMENTS  OR  ROTARY  INERTIAS  (REAL) 


TERMINATE  WITH  A ZERO  DATA  SET 


I 


-ELEMENT  LfTAD  MULTIPLIERS- 

ENTER  l DATA  SETS  < ELM ( I ># I ■ l * 4 ) WHERE  ELM( I ) ARE  THE  LOAD 
FACTORS  FOR  LOAD  CASES  (A.B.C.D) 

Q..Q..O..O. 


DO  YOU  WANT  A SUMMARY  OF  THE  SAP  INPUT  DATA  (YES  OR  NO  ) 
YES 


****  INPUT  DATA  SUMMARY  FOR  S.A.P.  *♦** 


NIFORMLY  LOADED  CANTILEVERED  REAM 


MASTER  CONTROL  INFORMATION- 


NUMNP  = 7 

NELTYP=  1 
LL  = l 

NF  = 0 

NDYN  = 0 

MODEX  = 0 

NAD  = 0 

KEQB  = 0 


-NODE  POINT  DATA- 


NODE 

NO 

COORD 

ID 

XI 

COORDINATES 

X2 

X3 

NODE 

INC 

NODAL 

TEMPERATURE 

1 

.0000 

.0000 

.0000 

0 

.0000 

2 

• 2000+0 1 

.0000 

.0000 

0 

.0000 

6 

.1000+02 

.0000 

.0000 

1 

.0000 

7 

.0000 

.1000+01 

• 0000 

0 

.0000 

NODE 

PERMANENT 

SINGLE 

POINT 

CONSTRAINTS 

NO 

XI 

X2 

X3 

R1 

R2 

R3 

1 

1 

1 

1 

1 

1 

1 

2 

1 

0 

1 

1 

I 

0 

6 

l 

0 

1 

1 

1 

0 

7 

1 

I 

1 

1 

I 

1 

-ELEMENT  DATA 
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THRKE-DIMENS I ONAL  BEAM  ELEMENT- 


NELMT- 

5 

NPROP= 

1 

NFEF  = 

0 

NMAT  * 

1 

MID 

L 

NU 

RHO 

WT 

1 

. 1004-08 

. '100 

.000 

.000 

PID 

AHEA 

SHEAR  1 

SHEAR  8 

J 

1-22 

1-33 

1 

t .00 

.000 

.000 

• 

O 

o 

1 .00 

1 .00 

GRAY 

ELEMENT 

LnAD  FACTORS 

DIR 

CASE  A 

CASE  B 

CASE  C 

CASE  D 

X 

.000 

.000 

.000 

.000 

Y 

.000 

.000 

.000 

.000 

Z 

.000 

.000 

.000 

.000 

BEAM 

ELMT 

DATA 

NODE 

NODE 

NODE 

MID 

PID  FIXED- END 

FORCE 

I.D. 

CODE  CODE 

GEN 

NO 

I 

J 

K 

A R 

C 

D 

I J 

PARAM 

1 

1 

2 

7 

1 

1 0 0 

0 

0 

0 0 

1 

5 

5 

6 

7 

1 

1 0 0 

0 

0 

0 0 

0 

- CONCENTRATED  L C1ADS /MASSES- 


NUDE 

NO 

LOAD 

CASE 

X 

FORCES 

Y 

Z 

X 

MOMENTS 

Y 

z 

2 

1 

• 000 

.2004-00 

.000 

.000 

.000 

.000 

3 

1 

.000 

.2004-00 

.000 

.000 

.000 

.000 

4 

t 

.000 

.2004-00 

.000 

.000 

.000 

.000 

5 

1 

.000 

.2004-00 

.000 

.000 

.000 

.000 

6 

I 

.000 

.2004-00 

.000 

.000 

.000 

• 000 

-ELEMENT  LOAD  MULTI PL  I ERS- 

LOAD  ELEMENT  MULTIPLIER 

NO  CASE  A CASE  B CASE  C CASE  D 

1 .000  .000  .000  .000 


IS  THERE  AN  ADDITION  CASE  (YES  OR  NO) 
NO 


NORMAL  EXIT.  EXECUTION  TIME: 


3258  MILLISECONDS. 
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4.2. 4. 2 Data  File  Mode 

The  input  data  for  the  data  file  mode  is  to  be  stored  in  an  element  called 
"EXAMPLE1".  The  editing  session  used  to  create  this  file  is  as  follows: 


!6D2§SA8§/Ti-U  :29-  (00 ) t F 
INPUT 

CANTILEVERED  BEAM  WITH  UNIFORM  LOADING 
7» 1 # 1 #0#0 

I>lll<t»l>l>lll*0a>0>>0w0l0> 

2,5.1 
1 *0 

1 / 1 • E7  * • 3#  0 • # 0 • 

1 ,1 .«0.»0«»1  ..i .» 1 • 

YES 

I#I>2»7»1»1»0»0»0#0»0»0»1 
1 #7 

5»5>6#7  » 1 » 1 »0»0»0»0»0>0»0 

2 , \ , 0 • , •2f0*>0*^0**0» 

3 # 1 >0>i  •2»0«*0»#0»#0« 

A , 1 , 0 • , •2*0.»0»>0«»0. 

5/l*0«»»2>0«»0»»0*»0» 

6#  1 ,0  9, 

0>0<0>i0o0*<0*i0*<0> 

o.»o.»o.»o. 

YES 

NO 
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After  verifying  that  the  data  file  is  correct.  The  I DP  is  executed  using 


the  element  EXAMPIE1  as  follows: 

• XQT  3YMP0S IUM*PREPROCESS OR • SAP 


*****  TS0006**  **  + 

INPUT  DATA  PREPROCESSOR  FOR  S.A.P. 

UNIVAC  1100-SERIES  VERSION.  INPUT  DATA  IS  TO  BE  ENTERED  FROM  EITHER 
THE  INTERACTIVE  OR  THE  DATA  FILE  MODE  USING  FREE  FORM  INPUT  AS 
DESCRIBED  IN  THE  UNIVAC  FORTRAN  V MANUAL.  IN  THIS  MODE  THE 
INDIVIDUAL  ELEMENTS  OF  THE  DATA  SET  ARE  SEPARATED  BY  COMMAS  (.). 
NOTE  THAT  THE  DECIMAL  POINT  (.)  IS  REQUIRED  FOR  REAL  NUMBERS • THE 
NORMAL  FORTRAN  CONVENTION  IS  FOLLOWED  FOR  VARIABLE  NAMES  UNLESS 
OTHERWISE  INDICATED  BY  APPENDING  AN  ANNOTATION  TO  THE  VARIABLE 
DESCRIPTION 


ENTER  NAME  OF  OUTPUT  DATA  FILE 
SAPIN 


DO  YOU  WANT  TO  ENTER  DATA  FRCM  A DATA  FILE  (YES  OR  NO) 
YES 


ENTER  NAME  OF  INPUT  DATA  FILE 
EXAMPLEl 


****  INPUT  DATA  SUMMARY  FOR  S.A.P.  **** 


CANTILEVERED  BEAM  WITH  UNIFORM  LOADING 


-MASTER  CONTROL  INFORMAT  I ON- 

NUMNP  = 7 
NELTYP"  1 
LL  * 1 
NF  = 0 
NDYN  a 0 
MODEX  ■ 0 
NAD  » 0 
KEQB  a 0 
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-NODE  POINT  DATA- 


NODE 

COORD 

COORDINATES 

NODE 

NODAL 

NO 

ID 

XI 

X2 

X3 

INC 

TEMPERATURE 

1 

.0000 

.0000 

.0000 

0 

.0000 

2 

.2000+00 

.0000 

.0000 

0 

.0000 

6 

.1000+02 

.0000 

.0000 

1 

.0000 

7 

.0000 

.1000+01 

.0000 

0 

.0000 

NODE 

PERMANENT 

SINGLE 

POINT 

CONSTRAINTS 

NO 

XI 

X2 

X3 

R 1 

R2 

R3 

1 

l 

1 

1 

1 

1 

1 

2 

1 

0 

1 

1 

1 

0 

6 

1 

0 

1 

1 

1 

0 

7 

1 

1 

1 

1 

1 

1 

-ELEMENT  DATA- 


- THREE-DIMENSIONAL  BEAM  ELEMENT 

NELMT-  5 
NPROP*  1 
NFEF  » 0 

NMAT  * I 


MID 

E 

NU 

RHO 

V T 

1 

•100+0H  .300 

.000 

.000 

PID 

AREA 

SHEAR 

1 

SHEAR  2 

J 

1-22 

1-33 

1 

1 .00 

.000 

.000 

1 .00 

1.00 

1 .00 

GRAV 

ELEMENT 

LOAD  FACTORS 

DIR 

CASE  A 

CASE 

B 

CASE  C 

CASE  D 

X 

.000 

.000 

.on^ 

.000 

Y 

.000 

.000 

• 

.000 

Z 

.000 

.000 

• 

.000 

BEAM 

DATA 

ELMT 

NODE  NODE  NODE 

MID 

PID  FIXED- 

END  FORCE 

I.D. 

CODE 

CODE 

GEN 

NO 

I J 

K 

A 

B C 

D 

I 

J 

PARAM 

1 

1 2 

7 

1 

1 0 

0 0 

0 

0 

0 

1 

5 

5 6 

7 

1 

1 0 

0 0 

0 

0 

0 

0 
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-CONCENTRATED  LOADS/MASSES- 


NODE 

NO 

LOAD 

CASE 

X 

FORCES 

Y 

Z 

X 

MOMENTS 

Y 

7 

2 

1 

.000 

.200+00 

.000 

.000 

.000 

.000 

3 

I 

.000 

.200+00 

.000 

.000 

• 000 

.000 

A 

1 

.000 

.200+00 

.000 

.000 

.000 

.000 

5 

1 

.000 

.200+00 

.000 

.000 

.000 

.0  00 

6 

1 

.000 

.100+00 

.000 

.000 

.000 

.000 

NO 

-ELEMENT  LOAD  MULTI  PL  I ERS- 


LOAD 

ELEMENT 

MULTIPLIER 

NO 

CASE  A CASE  B 

CASE  C 

CASE  D 

1 

.000  .000 

• 000 

.000 

NORMAL  EXIT.  EXECUTION  TIME* 
DATA  IGNORED  - IN  CONTROL  MODE 

3397  MILLISECONDS 
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4.2.4. 3.  Formatted  SAP  Input  File 

The  formatted  SAP  input  file  may  be  listed  by  means  of  the  text  editor.  The 
session  required  to  list  the  output  file,  SAPIN  is  shown  below. 


@ ED  SAPIN. 

ED  28C  03/14-11 S15-(00>:F 
EDIT 
♦ P 50 

UNIFORMLY  LOADED  CANT  I LEVERED  BEAM 


7 

1 

1 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

.000 

.000 

.000 

.000 

2 

I 

0 

1 

1 

1 

0 

2.00 

.000 

.000 

.000 

6 

1 

0 

1 

1 

1 

0 

10.0 

.000 

.000 

.000 

7 

1 

1 

1 

1 

1 

1 

.000 

1.00 

.000 

.000 

2 

5 

1 

0 

1 

1 

.100+08 

.300 

.000 

.000 

1 

1 .00 

.000 

.000 

1.00 

o 

o 

• 

1.00 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

1 

l 

2 

7 

1 

1 

0 

0 

0 

0 

0 

0 

5 

5 

6 

7 

1 

1 

0 

0 

0 

0 

0 

0 

2 

1 

.000 

.200 

.000 

.000 

.000 

.000 

3 

l 

.000 

.200 

.000 

.000 

.000 

300 

4 

1 

.000 

• 200 

.000 

.000 

.000 

• 000 

5 

1 

.000 

.200 

.000 

.000 

.000 

.000 

6 

l 

.000 

.200 

.000 

.000 

.000 

.000 

0 

0 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

.000 

EOF  AT  LINE 

21 

*G  19 

6 

1 

.000 

.200 

.000 

.000 

.000 

.000 

*C  /.2/. 

1 / 

6 

1 

.000 

.100 

.000 

.000 

.000 

• 000 

*E 

F.ND  EDIT 

21 

LINES 

OUTPUT 
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4.3.  SAP-Analysis  Program 

The  SAP  analysis  program  can  be  executed  in  either  of  two  modes,  both  of 
which  utilize  the  formatted  input  data  file  produced  by  the  IDP.  The  SAP 
analysis  program  can  be  executed  in  the  "Demand"  mode  for  small -to- intermediate 
sized  problems;  and,  for  large  problems  may  be  submitted  to  the  "Batch"  queue. 
The  only  important  difference  in  the  two  modes  is  that  the  user  is  free  to 
perform  other  tasks  on  the  terminal  if  the  batch  mode  is  utilized. 

4.3.1.  Demand  Mode 

The  SAP  analysis  program  is  executed  by  entering  the  following  system 
commands  from  the  TTY; 

,8XQT  SYHPQS I UM*PROGRAM. SAPI 1 1 
@AUD  <SPEC> 


where  <SPEC>  is  the  name  of  the  input  file  which  was  created  by  the  preprocessor. 
4.3.2.  Batch  Mode 

The  SAP  program  may  be  submitted  to  the  "batch"  queue  by  the  command 
(3BATCH  SYMPYPOSIUM*PROGRAM.SAPIII,  <data>  , RMENTR 

where 


<data>  is  the  name  of  the  data  file 
created  with  the  preprocessor. 

4.3.3.  Redirecting  Output 

You  may  choose  to  have  tne  program  output  saved  on  a cataloged  file  in 
either  the  demand  or  batch  mode.  In  this  case  you  must 

(1)  Catalog  a file 

(2)  Inform  the  system  to  redirect  output 

(3)  Execute  the  program 

(4)  If  on  demand,  inform  system  to  direct  output  to  TTY  after 
execution 
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The  basic  commands  necessary  to  perform  steps  (1),  (2)  and  (4)  are  given 
in  section  2. 

The  third  step  depends  on  the  mode  of  operation.  The  following  are  examples 
of  run  stream  controls  for  both  modes. 

4. 3. 3.1.  Demand 

SCAT,  U SAPOUT 
0FREE  SAPOUT. 

3BRKPT  PRINT$/SAPOUT 

tfXQT  SYMPOS IUM*PROGRAM. SAP 1 1 1 

0ADD  SAP Id. 

OBRKPT  PRIWT$ 

At  this  point  control  would  be  returned  to  the  remote  terminal,  and  the 
SAP-output  could  be  examined  using  the  text  editor. 

4. 3. 3.2  Batch 

In  order  to  redirect  the  output  to  a cataloged  file  in  the  batch  mode  we 
must  create  a procedure  file  and  use  the  system  command  0START  rather  than  0BATCH 
to  submit  the  job  to  tne  batch  queue  as  indicated  in  section  2. 3.8.2.  The 
procedure  file  is  created  using  the  text  editor.  In  the  example  shown  below  the 
file  START  is  created  which  contains  the  required  run  stream  card-images. 

These  images  are  then  submitted  to  the  batch  queue  by  the  command  OSTART  START. 
3ASG.U  START, F (assign  file  "START") 

@EU  START.  (use  editor  to  create  file) 

#RUN  <RYBUD> , <ACC0UNT  N0>,  <PR0JECT>,  <TIME>,  <PAGES> 

#CAT,U  SAPOUT, F 
#FREE  SAPOUT. 

IBRKPT  PRINTS/ SAPOUT 


4-17 


#AL)0  SAPIN. 

#FIM 

C-R  (Return  to  edit  mode) 

#C,*  /#/3/  (change  #'s  to  3's) 

3ST ART  START. 

dote  that  3-signs  cannot  be  entered  directly  in  the  editor  input  mode. 
The  *'s  have  thus  been  entered  in  place  of  the  3-sign  and  the  editor  nas  jeen 
utilized  to  change  # to  3. 

4.3.4.  Scanning  the  output 

The  text  editor  may  be  used  to  scan  the  output  by  entering 
3ED  SAPOUT. 

The  displacements  may  be  found  by 
L DIST 

The  forces  may  be  found  by 
L FORCE 

The  lines  of  text  following  the  selected  current  line  may  be  printed  by 
0 ML 

Where  NL  is  the  number  of  lines  to  be  printed. 

4.3.5.  Printing  Output 

The  output  may  be  printed  on  the  high-speed  printer  by  entering 
3SYM  SAPOUT.,, RMENG 

If  you  do  so,  and  wish  to  have  the  output  mailed  you  must  use  the  text 
editor  to  enter  your  name  and  mailing  address  in  the  output  file. 


L 
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4.3.6.  Analysis  Program  Output. 


The  output  of  the  SAP-analysis  program  for  the  example  problem  is  as 
follows. 
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Table  4.1-SAP  Preprocessor  Input  Data  Elements 


1 


I.  72  Column  Title 

II.  (NUMNP .NELTYP ,LL ,M0DEX ,KEQB ) 

If  LL.LE.O  Go  to  IIA 

GT.O  Go  to  III 
IIA.  (NF,KDYN,ND) 

III.  Data  Sets 

(IC,N,IX(1  6) ,(X,Y ,Z) ,KN,T) 

MUST  END  WITH  N= NUMNP 

IV.  (NTYPE ,NELMT ,NMAT ) 

A NON-ZERO  ENTRY  FOR  NMAT  REQUIRED  FOR  ALL  ELMT.  TYPES 
IV. 1 FOR  NTYPE. EQ.l  IV. 1 

.A  NMAT  SETS  (MID.E.a.p.A.y) 

.B  SKIP  IF  LL.EQ.O 

YES  IF  ALL  ELMT  LOADS  ARE  ZERO 
4 SETS  (A,B,C,D)  OTHERWISE 
.C  DATA  SETS 

(NEL,I,J,MID,T,K) 

MUST  END  WITH  NELCNELMT 
IV. 2 FOR  NTYPE. EQ. 2 
.A  (NPROP.NFIX) 

.B  NMAT  SETS  (MID.E.v.p.y) 

.C  NPROP  SETS  (PID,A,SA2,SA3,J,I22,I33) 

•D  SKIP  IF  LL.EQ.O 

YES  IF  ELMT  LOADS  ARE  ZERO 
3 SETS  (A,B,C,D)OTHERWISE 


L 
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TABLE  4-1  (continued) 


.E  SKIP  IF  NFIX.EQ.0 


NFIX  SETS  (NFIX,FEFATI(1  6),FEFATJ(1  6)) 

.F  DATA  SETS  (NELM,I,J,K,MID,PID,IDA,IOB,IOC,IDD,IREL,JREL,KG) 
MUST  END  WITH  NELM=NELMT 
IV. 3 FOR  NTYPE.EQ.3 
.A  (NTEMP, MODE) 

.B  NMAT  SETS  AS  FOLLOWS: 

,B .1  (MID, NTEMP .y.p.B) 

.B.2  NTEMP  DATA  SETS  (T ,En ,ES  ,Et  ,x>ns  ,vnt ,v$t  ,Gn$  ,<*n ,a$ ,«t) 
.C  SKIP  IF  LL.EQ.O 

YES  IF  ELMT  LOADS  ARE  ZERO 
4 SETS  (TL.PL.XL.YL.ZL)  OTHERWISE 
.D  DATA  SETS  (NEL ,1 ,J ,K,L .MID.TREF.PRES .OP.KG.T) 

MUST  END  WITH  NEL»NELMT 
IV. 4 FOR  NTYPE.EQ.4 

.A  ( NTEMP, NOPT .MODE) 

.B.l  (MID,NTEM,y,p,B) 

.B.2  NTEM  SETS  (T.En,Es,Et,vns,vnt,vst,Gns,an,as,at) 

.C  SKIP  IF  LL.EQ.O 

YES  IF  ELMT  LOADS  ARE  ZERO 
4 SETS  (TL.PL.XL.YL.ZL)  OTHERWISE 

• D DATA  SETS  (NEL, I .J.K.L.MID.TREF.PRES.NOP.KG.T) 

WST  END  WITH  NEL*NELMT 
IV. 5 FOR  NTYPE.EQ.5 


.B  NMAT  SETS  (MID.E.v.Y.a) 
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TABLE  4-1  (continued) 


.C  SKIP  IF  LS.EQ.O 

LS  SETS  (LID.LT.P.Y.NFACE) 

.D  G 

.E  SKIP  IF  LL.EQ.O 

YES  IF  ELMT  LOADS  ARE  ZERO 
5 DATA  SETS  (A.B.C.D)  OTHERWISE 
.F  DATA  SET  (EID,N0DESO-8),I0,MID,INC,LS(1-4)  ,FACE(1  ,2)  ,T) 
MUST  END  WITH  EID=NELMT 
IV. 6 FOR  NTYPE.EQ.6 

.A  ANSWER  TO  "IS  MATERIAL  ISOTROPIC?" 

.A.l  YES,  NMAT  SETS  (MID.P.a.E.v)  .0^1-5 
.A. 2 NO,  NMAT  SETS  (MID,P,«x,ay ,«z,Cxx’Cxy ,CxS ,Cyy,Cys ,Cxy) 
.B  SKIP  IF  LL.EQ.O 

YES  IF  ELMT  LOADS  ARE  ZERO 
5 DATA  SETS  (A.B.C.D)  OTHERWISE 
.C  DATA  SETS  (NEL.I .J.K.L.NO.MID, KG ,T  .PRESS .TEMP .GRAD) 

MUST  END  WITH  NEL*NELMT 
IV. 7 FOR  NTYPE.EQ.7 

.A  SKIP  IF  LL.EQ.O 

YES  IF  ELMT  LOADS  ARE  ZERO 
DATA  SET  (A.B.C.D)  OTHERWISE 
.B  DATA  SETS  (NODE.I.J.K.L.ID.IR.KN.DISP.ROT.STIFF) 

LAST  ELMT  MUST  BE  DEFINED 
IV. 8 FOR  NTYPE.EQ.8 
.A  NL 

.B  NMAT  SETS  (MID,E,NU,Y,«) 

.C  SKIP  IF  NL.EQ.O 

NL  DATA  SETS  (LID.LT.P.Y.NFACE) 
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TABLE  4-1  (continued) 


I 

I 

1 

I 


L 


.D  (TREF.G) 

.E  SKIP  IF  LL.EQ.O 

YES  IF  ELMT  LOADS  ARE  ZERO 
5 DATA  SETS  (A,B,C,D)  OTHERWISE 
.F  DATA  SETS  (EID,I0,MID,INC,LS(l-*4)  ,T,H0DE(W16)} 

MUST  END  WITH  EID=NELMT 

IV.  9 FOR  NTYPE.EQ.9 

.A  (MAXTP.NSECT.NBRP.MAXTAN) 

.B  NMAT  SETS  AS  FOLLOWS: 

.B.l(M.NT) 

.B.2  (LABEL (36  cols.)) 

.B.3  NT  DATA  SETS  (T.E.v.a) 

.C  NSECT  SETS  AS  FOLLOWS: 

.C.l  (N,D,T ,SF,WT,RHO) 

.C.2  (LABEL (18  cols.)) 

.D  NBRP  NODE  NUMBERS 
.E  SKIP  IF  LL.EQ.O 

YES  IF  ELMT  LOADS  ARE  ZERO 
5 DATA  SETS  (A.B.C.D)  OTHERWISE 
.F.l  DATASETS  (N.CODE.I  .J.MID.SID.T.P^  ,A2,A3,K) 

.F.2  OMIT  IF  CODE.EQ.O  (R,CODE,X,Y,Z,F) 

MUST  END  WITH  N»NELMT 

V.  (N,L,FX,FY,FZ,MX,MY,MZ) 

TERMINATE  WITH  ZERO  SET  1 

VI.  SKIP  IF  LL.EQ.O 

LL  DATA  SETS  (A.B.C.D) 

IF  (LL.GT.O)  ANSWERS  TO  FOLLOWING  QUESTIONS,  OTHERWISE  PROCEED  TO  VII. 
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TABLE  4-1  (continued) 


Q.l  “DO  YOU  WANT  A SUMMARY  OF  SAP  INPUT  DATA?"  (YES  OR  NO) 
Q.2  "IS  THERE  AN  ADDITIONAL  CASE?"  (YES  OR  NO) 


VII.  (IFNDYN.EQ.l  ,2,  or  3)  (IFSS.NITEM.RTOL.COFQ.NFO) 

.A  (NDYN.EQ.2  or  4) 

.A.,  (NFN 'NGM.NAT ,HT .NOT ,OT , * ALPHA. BETA  <!f  (SSJn ! lo ! 4 )> 
.A. 2 SKIP  IF  (NFN.EQ.O)  (N.IC.IFN.IAT.XM) 

TERMINATE  WITH  A ZERO  DATA  SET 
.A. 3 SKIP  IF  NGM.EQ.O  (NFNX,NFN,Y,NFN,Z,NATX,NAT.Y .NAT.Z) 

.A. 4 SKIP  IF  NAT.EQ.O 

NAT  ARRIVAL  TIMES  (A(I) ,1*1  ,NAT) 

.A. 5. A NFN  SETS 
(NLP.SFTR) 

HED  (LE.60  COLS) 

.A.5.B  NLP  PAIRS  (T,F) 

.A.6.A.1  (KKK.ISP) 

.A.6.A.2  (N,IC(I),I*1  ,6) 

TERMINATE  WITH  ZERO  SET 
.A.6.B.1  (KKK.ISP) 

.A.6.B.2  (N,IS(I) ,1*1 ,12) 

TERMINATE  WITH  A ZERO  SET 

Answer  Q.l  and  Q.2,  above,  to  terminate  dynamic  input  data 
.B  (NDYN.EQ.3) 
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5.0  VISCEL/ELAS  - Viscoelastic  Analysis  of  Linear  Systems 

5.1  Introduction 


VISCEL,  a general-purpose  computer  program  for  the  analysis  of  linear 
viscoelastic  structures  is  described  in  reference  (5).  This  program  contains 
the  ELAS-program  capability  as  a subset;  and  incorporates  the  basic  ELAS  input 
data  set.  The  ELAS-program  and  input  data  are  described  by  reference  (6).  The 
modifications  to  the  basic  ELAS  input  required  to  perform  a viscoelastic  analysis 
are  described  in  reference  (5). 

The  VISCEL  and  ELAS  programs,  together  with  program  documentation,  are 
available  from  COSMIC  (7).  The  ELAS  program  is  also  maintained  and  distributed 
by  its  author;  and,  is  available  from  the  Duke  University. 

This  report  describes  the  execution  of  the  (1)  an  interactive  program  which 
generates  the  formatted  input  file  for  VISCEL;  and  (2)  the  execution  of  the  input 
data  processor  (IOP)  and  the  VISCEL  analysis  program  on  the  University  of 
Maryland  Uni  vac  1103. 

5.2  VISCEL/ELAS  Input  Data  Preprocessor 

The  VISCEL  IDP  operates  in  the  Demand  mode  and  allows  the  user  to  input 
unformated  data  in  a program-control  led  logical  sequence.  The  IDP  performs 
some  elementary  error  checks  and  writes  a file  which  contains  the  formated  input 
required  by  the  VISCEL  analysis  program.  The  IDP  operates  in  either  of  two 
modes,  interactive  or  data  file;  and,  further  allows  the  user  to  control  the 
amount  of  computer  generated  prompting  in  the  interactive  mode. 

, In  the  interactive  mode  the  IDP  Issues  prompting  information  which  is 

displayed  by  the  terminal  device.  The  user  then  inputs  the  requested  informa- 
t tion  directly  by  means  of  the  keyboard  in  unformated  form;  i.e.  the  data  elements 

I 
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are  delimited  by  commas.  This  interaction  continues  until  all  of  the  input 
requirements  have  been  satisfied.  In  the  data  file  mode  the  user  prepares  the 
unformatted  input  to  the  IDP;  and,  then  informs  the  IDP  the  name  of  the  data 
file.  The  IDP  then  reads  the  input  data  directly  from  the  data  file  in  exactly 
the  same  sequence  as  the  interactive  mode;  but,  completely  bypasses  the  prompting 
information. 

In  either  mode  the  sequence  required  for  the  IDP  input  is  shown  by  Table  5.1. 
For  tne  most  part  the  input  data  elements  have  the  same  names  and  meaning  as 
those  variables  defined  in  the  program  users  manual.  However  it  is  recommended 
that  the  interactive  mode  be  executed  initially  with  maximum  promoting  in  order  to 
completely  define  each  data  element. 

The  user  will  note  that  the  IDP  does  not  request  all  of  the  program  control 
information  described  in  the  ELAS/VISCEL  users  manuals.  The  IDP  does  not  allow 
the  use  of  externally  prepared  subroutines  such  as  "CORG",  "MESG",  and  "BUNG", 
nor  does  it  allow  the  inclusion  of  data  cards  such  as  those  required  by  the 
relaoelling  options,  ISHUF  = 2 or  3.  The  following  parameters  are  preset  in 
the  IDP. 

ITAP  = 2 
ITAS  = 9 
ICOR  = 0 
I BUN  * 0 
IMES  * 0 

In  addition  the  following  parameters  have  a more  restrictive  definition 
than  given  in  tne  users  manual: 

0 < ISHUF  < 1 
0 < IPIR  < 1 
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The  user  will  notice  that  the  IDP  requests  redundent  information.  For  example, 
the  parameter  IT  is  requested  as  control  information;  and,  the  number  of  elements 


I 

I 

is  also  required  in  element  descriptions.  This  redundancy  is  required  in  order 
to  simplify  tne  IDP;  if  tnese  redundancies  were  eliminated  the  data  management 

I 

tasks  would  be  increased.  The  IDP  does  contain  an  error  check  to  assure  that 

redundant  control  parameters  are  consistant. 

You  will  note  that  the  IDP  requires  a data  item  in  addition  to  that  normally 

required  by  the  analysis  program  - the  class  of  structure.  The  various  structural 

classes  are  described  by  the  interactive  mode  of  the  IDP  with  maximum  prompting; 

and  are  defined  by  Table  A-5  of  reference  (5).  The  class  number  determines  tne 

type  of  elements  that  may  be  used;  and,  the  IDP  logic  restricts  the  user  to  the 

class/element  relation  defined  by  the  table  referenced  above. 

Tne  only  difference  between  a normal  ELAS  run  and  a viscoelastic  run  is  the 

inclusion  of  the  parameter  ISUCA  which,  again  is  defined  in  the  IDP  as  follows: 

. EQ.O  Signals  termination  of  the  input;  if  ISUCA  has  not 
been  previously  defined  then  this  indicates  that  a 
standard  ELAS  run  is  to  be  made 

. LT.O  Allows  multiple  runs 
. GT.O  Indicates  Viscoelastic  run 
. EQ . 1 Master  data  and  following  modified  data 
. GT.l  For  additional  modified  data  sets. 
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5.2.1  Interactive  Mode 

In  the  Interactive  mode  the  IDP  sets  up  a question  - answer  dialogue  with 
the  user.  The  IDP  defines  a data  set  to  be  Input  and  then  accepts  the  Input  In 
unformatted  form  directly  from  the  terminal  keyboard.  The  ELAS/VISCEL  IDP  has 
two  level  of  prompting  which  Is  controlled  by  user  request.  The  minimum 
prompting  mode  requests  data  elements  by  displaying  the  data  element  mneumonics; 
while  the  maximum  prompting  mode  defines  each  of  the  data  elements. 

In  either  the  Interactive  or  the  data  file  mode  the  IDP  creates  a formatted 
Input  file  for  the  VISCEL/ELAS  analysis  program.  This  file  which  is  stored  on 
a mass  storage  device  must  be  identified  by  means  of  a file  name.  The  program 
thus  requests  the  user  to  enter  a name  for  the  file.  The  IDP  will  then  catalog 
the  file;  and,  an  error  will  result  if  the  file  already  exists.  The  user  should 
note  that  the  file  is  cataloged  as  a private  file  and  may  be  named  using  read 
keys  as  follows. 

QUAL I F I E RENAME/ KEY 1 / KE  Y2/ 

where  KEY1  is  the  write  key;  KEY2  is  the  read  key.  The  user  is  cautioned  to 
retain  these  keys  in  a secure  place,  if  used. 

5.2.2  Data  File  Mode 

In  this  mode  the  user  prepares  a data  file  containing  the  responses  to  the 
IDP  either  off-line  using  a paper  tape  or  tape  cassette;  or  on  line  using  the 
text  editor.  The  data  elements  are  described  by  Table  5.1.  The  purpose  of  this 
input  mode  is  to  allow  the  user  to  prepare  an  input  file  for  the  IDP  so  that  the 
prompting  produced  by  the  IDP  can  be  bypassed. 

The  IDP  will  query  the  user  to  determine  if  the  data  file  mode  is  to  be 
used.  If  so,  the  user  is  requested  to  enter  the  name  of  the  file  containing  the 
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prepared  file  of  unformatted  input.  The  IQP  will  then  read  from  this  file  to 
satisfy  read  requests  completely  by-passing  the  interactive  prompting. 

Occasionally  the  user  may  make  an  error  in  preparing  the  data  file  which 
will  cause  a read-error  in  the  IOP,  and  an  error  exit.  The  information  printed 
by  the  system  at  this  time  will  be  of  little  assistance  in  identifying  the  data 
elements  which  are  in  error.  A recommended  procedure  at  this  point  is  to  place 
the  IDP  in  the  prompting  mode;  but  to  cause  the  IDP  to  read  the  pre-created  data 
file.  This  can  be  accomplished  by  responding  "NO"  to  the  question  "00  YOU  WANT 
TO  ENTER  DATA  FROM  A DATA  FILE?".  Then  in  response  to  the  first  prompt  for  data 
the  data  file  is  added  by  the  command 

@ADD  QUALIFIER*NAME 

The  IDP  now  prints  all  prompting,  but  reads  the  data  from  the  data  file. 

It  is  then  possible  to  identify  the  data  elements  in  error  by  noting  at  which 
point  the  error  exists. 

5.2.3  IDP  Input  Requirements 

The  order  in  which  the  IDP  prompts  the  user  for  input  differs  slightly 
from  the  ELAS/VISCEL  manuals,  and  provides  default  values  for  scratch  tapes  etc. 
In  addition,  the  IDP  does  not  allow  the  user  to  input  parameters  which  would 
require  user  prepared  subroutines. 

5.3  Examples 

The  ELAS  examples  are  taken  from  reference  8. 

5.3.1  Planar  Truss  Case 

The  planar  truss  is  as  shown  by  figure  5.1 
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FIGURE  5.1  Planar  Truss 


Before  proceeding  to  execute  the  IDP  we  will  have  to  refer  to  the  ELAS/VISCEL 
manuals  to  define  certain  of  the  control  parameters.  Referring  to  Table  5.2, 
which  is  taken  from  reference  (5)  we  see  that  the  planar  truss  is  structural 
class  number  one;  and,  that  there  Is  only  one  element  type  associated  with  this 
structural  class;  the  planar  truss.  Finally,  referring  to  Table  5.3  we  see  that 
the  parameters  IDEG,  NVERT,  NWDS  have  the  following  values: 

IDEG  - 2 (column  4 of  Table  5.3) 

NVERT  = 2 (column  3 of  Table  5.3) 

NWDS  = 5 (column  5 of  Table  5.3) 

With  this  information  we  can  write  the  data  elements  of  data  sets  2,3,4,  and 
5 as  follows: 

Set  number  Data  Elements 

2 (11,20,2) 

3 (0,  0,  0,  2,  5,  4) 

4 (1,0,  1,  0,  0,  0,  0) 

5 (2,  0,  0,  0,  0) 


Where  we  have  indicated  by  setting  NDBC  = 4 that  four  displacement  boundary 
conditions  are  to  be  supplied;  and  by  ITYPE  * 0 that  isotropic  material  properties 
are  to  be  used. 
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5. 3. 1.1  Interactive  Mode 

The  interactive  input  mode  for  the  creation  of  the  formatted  ELAS/VISCEL 
input  file  is  as  follows;  where  the  formatted  ELAS  input  data  is  to  be  written 
onto  a file  called  PR0BO1 . 

SXOT  SYMPOSIUM*PREPROCESSOR.  VI  SCF1. 

*****  TS007  ***** 

INPUT  DATA  PREPROCESSOR  FOP  ELAS/VISCFU 

fc^JTFR  NAME  OF  riUTPUT  DATA  FILF 
PROPO  1 


Dr i YOU  WANT  TO  ENTER  DATA  FROM  A DATA  FILE  ( YFS  nR  MO) 

MO 


DO  YOU  WANT  MINIMUM  PROMPTING  (YFS  OR  NO) 
NO 


UNI  VAC  1100-SF.RIES  VERSION.  INPUT  DATA  TO  BE  ENTERED  FROM  FITHFr 
THE  INTERACTIVE  MODE  OR  THE  DATA  FILE  MODE  USING  FREE  FnRM  INPUT  AP 
DESCRIBED  IN  THE  UNIVAC  FORTRAN  V MANUAL.  IN  THIS  MODE  THF  INDIVIDUAL 
Fl_  FM  FN  T S OF  THE  INPUT  DATA  SET  ARE  SEPARATED  BY  COMMAS.  THE  DECIMAL 
POINT  (.)  IS  REQUIRFD  FOR  RFAL  NUMBERS)  THE  NORMAL  FORTRAN  SIGN  CONVFN 
SrON  IS  FOLLOWED  FOR  VARIABLE  NAMES  UNLESS  OTHERWISE  NOTED. 


ENTER  TITLE  (72  COLUMNS  OR  LESS) 
BLANF  TRUSS  PROBLEM  IN  MANUAL  Nn  8 


- CONTROL  INFORMAT I ON- 

ENTF:R  THF  DATA  SET  ( I N,  I T,  I DFG  ) 

IN  = TOTAL  NUMBER  OF  NODES 
IT  = TOTAL  NUMBER  OF  FLEMENTS 

I DEG  = NUMBER  OF  DEGREES  OF  FREEDOM  AT  NODE  POINTS 

1 1,20,  2 
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j 


ENTER  THE  DA  I A SET  C I TYPF  .,  I G FM,  I STR,  M VERT,  NWDS,  NDBC  ) 
I TYPE- MATER I AL  INDICATOR 
. EO.O  ISOTROPIC 
. EU.l  URTHI  'TROPI  C 
.EQ.2  ANISO TROPIC 
IC  EM  -GEOMETRY  INDICATOR 

.EO.O  TWO-1  IN  X- Y PLANE 
• K.O  • 1 THRE  E - D 

I^TR  -PLANE;  STRAIN  INDICATOR 
. EQ.l  PLANE  STRAIN 
.EO.O  PLANE  STRESS 
NUEin=MA*  NO  OF  ELEMENT  NODES 

NUDS  = MA'<  NO  OF  INPUT  ITEIMS  TO  DESCRIBE  ELEMENT 
NDBS  =TOrAL  NO  OF  DEFLECTION  BOUNDARY  CONDITIONS 

0»  0)  Q>  3>  5,  A 


ENTER  THE  DATA  SET 
IMAT  = TOTAL  NO 
NT I C = TOTAL  NO 
I AREA=  TOTAL  NO 
I TOP  = TOTAL  NO 
IYY  = TOTAL  NO 
IZZ  = TOTAL  NO 
IMFI  = TOTAL  NO 
2,  Q,  3/  0,  0 , 0,  U 


( IMAT,NTIC,  IARF.A,  I T I IR,  IYY,  IZZ,  IMFI  ) 

OF  MATERIALS 

OF  THICKNESSES 

OF  CROSS- SECTI  rNAL  AREAS 

OF  TORSIONAL  CONSTANTS 

OF  MOMENTS  OF  INERTIA  ABOUT  Y-AXIS 

OF  MOMENTS  OK  INERTIA  ABOUT  7-AXIS 

OF  ANGLES  FIXING  LOCAL  Y AND  7-AXFS 


ENTER  THE  DATA  SET  ( IP,  I PRS,  I TEMP,  I G RADY  , IG  RATT  ) 

IP  = T OT  AL  NO  OF  CONCENTRATED  LOADS 
IPRS  = TOTAL  NO  OF  PRESSURES 
ITEMP  = TOTAL  NO  OF  TEMP  INCRFIASES 
IGRADY=TOTAL  NO  OF  TEMP  GRADIENTS  ALONG  LOCAL  Y 
IGRAD7,=  T0TAL  NO  OF  TEMF  GRADIENTS  ALONG  LOCAL  Z 
8, 0, 0, 0, 0 


ENTER  LINK  AFTER  WHICH  EXECUTION  IS  TERMINATED 
.EQ.l  INPUT 
.EQ.2  GENERATION 
.EG.  3 DEFLECTION 
.EQ.4  STRESS 
A 


ENTER  PRINT  PARAMETER 
•EO.O  MINIMUM 
.EQ.l  INTERMEDIATE 
.EQ.2  DETAILED 


1 


ENTER  RELABEL  I MG  PARAMETER 
.EO.O  MO  RELABEL 
• EQ  • 1 RELABEL 


ENTER  COORD  SELECTION  INDICATOR  FOR  SH FLL S 

• EO.O  LOCAL  X IS  1-2  LIME  OF  LOWEST  NUMBERED  ELMT 
. EQ.l  PRINCIPAL 


ENTER  (G  1#  G2,G3)-DIRECTI  ON  COSINES  OF  ACCELERATION 

0 « j 0 • » 0 • 


ENTER  MAGNITUDE  OF  ACCELERATION  VECTOR  TIMES  IN  I T MASS 

0. 


-MATERIAL  PROPERTIES- 


ENTER  2 DATA  SETS  CE,G,  ALPHA) 

E =M ODLL US  OF  ELASTICITY 

G = SH  EAR  MODULUS 

AL PHA=  CUEF  OF  THERMAL  EXPANSION 
MATERIAL  NO  1 
2.50  OOP  » »0»>Q. 

MATERIAL  NO  2 
500000 .,0..0. 


-CROSS  SECTIONAL  AREAS- 

ENTER  3 DATA  SETS  Cl, AREA) 

I = AREA  SET  I.D. 

AREA= CROSS- SECTIONAL  AREA 

U-U 

THF  INTERPRETATION  OF  MEANINGLESS  INPUT  WAS  ATTEMPTED. 

THE  FOLLOWING  RECORD  IS  ERRONEOUS  OR  DOES  NOT  CORRESPOND  TO  FORMAT  SPF.CI 
FICATIONS: 

F,  .3 

I/O  CALLED  AT  SEQUENCE  NUMBER  000133  OF  TS010 


Notice,  that  the  ry&tem  hat  made  an  error  in  interpreting  the  input  data. 
Rather  than  terminating,  the  tertian  wilt,  be  continued ; and,  the  text 
editor  wilt  be  utilized  to  modify  the  formatted  data  bile  PR08OJ. 
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1 


-NODE  pqini  geomethy- 


IS  HIS  A PI. AMAH  PROBLEM  CYES  OH  NO) 
YES 


ENTER  11  DATA  SETS  <X,Y> 
modi-  1 

0  . . 0 . 

MODE  3 

fi..  in. 

NODE  3 

n . . g n . 

NODE  A 

1 0. , n «- 1 0. 

NODE  S 
1 0. < 20. 

NODE  6 

1 5. . 1 S. 

NODE  7 

2 0 • . 1 0 • 

NODE  H 

NODE  9 

NODE  10 
30 . . 2Q« 

NODE  1 1 

3Q*  r Qi 


-DEFLECTION  BOUNDARY  CnNDITIONS- 

ENTEH  ASETS  < NODE#  I DOF,  NODE-  P.  I DnF- P.  VALUE) 
MODE  =NODE  NO  OF  INDEPENDENT  D.  0.  F . 

I DOF  =D.O.F.  NO. 

NODE-P=NODE  NO  OF  RELATED  D.  O.F. 
IDOF-P*D.  O.F.  NO. 

VALUE  COEFFICIENT 

1.  1 # 1 . 1 . 0 » 

1 . 2.  1 . 2.  0« 

1 1. 1, 1 1.  l.Q. 

1 1.2.  1 1.2.0. 
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-ELEMENT  DESRIPTIONS 


THE.  TYPES  OF  FINITE  ELEMENTS  WHICH  'I  AY  BE:  USED  TO  PE'.SC'I  OF  ANY 
STRUCTURAL  SYSTEM  IS  DE'PENDE  \J  T ON  THE  CLASS  OF  STRUCTURAL  SYSTEM  • 

THE  CLASSES  HF  STRUCTURAL  SYSTEMS  ALL OWED  IN  THIS  PROGRAM  A'<F  AS 
FOLLOWS 

CLASS  DESRIPTIOM 

1 PLANAR  TRUSS 

2 SPACE  TRUSS 

3 PLANAR  FRAME 

A SPACE  FRAME 

5 GRIPWORK  FRAME 

6 PLANE  STRESS 

7 PLANE  STRAIN 

H PLATE  BENDING 

9 GENERAL  SOLID 

in  GENERAL  SHELL;  BENDING*  MEMBRANE 

11  GENERAL  SHELL*  MEMBRANE 

IS  SOLID  OF  REVOLUTION 

13  SHELL  OF  REVOLUTION*  MEMBRANE 

1 A SHELL  OF  REVOLUTION?  BENDING*  MEMBRANE 


ENTER  CLASS  OF  STRUCTURE 

1 


YOU  HAVE  SPECIFIED  STRUCTURAL  CLASS  1 

THE  COMPATABLE  ELEMENTS  FOR  THIS  STRUCTURAL  CL  ASA  A»-'E  ^PFCIFIEP  BY 
TABLE.  I I 1-2  OF  THE  ELAS  MANUAL?  AND*  YOU  ARE  RESTRICTED  TO  THESE 
ELEMENTS  BY  THE  PREPROCESSOR. 


ENTER  NUMBER  OF  TYPE  1 ELEMENTS 

20 


ENTER  20  DATA  SETS  OF  THE  E ALLOWING  FORM 
CM  1 * N2»  NMAT* I AREA* I PRESS* ITEMP) 

N 1 = NODE  NUMBER  OF  FIRST  VERTEX 

NS  =MODE  NUMBER  OF  SECOND  VERTEX 

MAT  =MATERIAL  TYPE  NUMBER 

NAREA  aCROSS-SE CTI ONAL  AREA  TYPE  NUMBER 

I PRESS=PRESSUE  TYPE  NUMBER 

HEMP  = UN  I FORM  TEMP  INCREASE  TYPE.  NUMBER 


ELMT 

i 

I 

*2*1* 

1*0* 

0 

? 

©ELMT 

2 

2 

*3*1* 

1*0* 

0 

ELM  T 

3 

1 

» 4#  1 » 

1*0* 

0 

\ 
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ELMT  A 

2,  A,  1/  l»0,0 
ELMT  S 

3, 4,  1 , 1 , Q, 0 
EL  -IT  6 

3,  5,  1,  2,0,  C 
ELMT  7 

A,  5,  1 , P,  0,  <• 

ELMT  « 

^)7< 1 , 2 » 0 , 0 

ELMT  9 

A,  6,  1 , 2, 0, 0 

ELMT  10 

2,  6,  1 , p>  0/  0 

ELMT  1 1 

ELMT  12 

6.7.2,  1,0,0 
ELMT  13 

S^Bj^UOiO. 
FLMT  14 

7.8.2,  1,0,0 
ELMT  15 

7.11,a.li0,'i 
ELMT  16 
7, 9,2, 3, 0,0 
ELMT  17 

Ltia*  2jJU^jJL 

ELMT  IB 

B,  1 0,  2,  3,  0,  0 
ELMT  19 

9,  1 1,2,  3,  0,  0 
ELMT  20 
9,  10,2,  3,  0,  0 


- CONCENTRATED  NODE  FORCFS  AMD  MOMF.NTS- 

EMT E.r<  2 DATA  SETS  OF  THE  FORM  (MODE,  I DOF,  VALUE)  WHERE: 
NODE-MODE  LABEL  MO. 

I DOF=  DEC  REE  OF  FREEDOM  CODE 
. EQ.l  FORCE  IN  X-DIR 

• EQ • 2 FORCE  IN  Y-DIR 

• EO . 3 FORCE  IN  Z-DIH 
.EQ.4  MOMENT  ABOUT  *-AXIS 
.F,Q.S  MOMENT  ABOUT  Y-AXIS 

• EQ*  6 MOMENT  ABOUT  7-AXIS 
VALUE-MAGM I TUDE  OF  FORCE  COMPONENT 

5.2, - 10* 

8. 2,  - 1 0.- 


ENTF.R  INCREMENTAL  PARAMETER  ISUCA 

. EQ.O  STANDARD  LINEAR  ELASTIC  (ELAS)  SOLUTION 
*LT.  1 MULTIPLE  RUN 

*GE.  1 VALUE  OF  M FOR  VISCOELASTIC  SOLUTION 
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DO  you  want  a summary  of  thk  input  data  <yfs  or  no> 

rFS 


summary  of  ELAS/VISCFL  INPUT  rata- 


plan k TRUSS  problem  in  manual  no  h 


-CONTROL  INFORMATION- 


NUMBER  OF  NODES 11 

NUMBER  OF  ELM TS 20 

D.O.F  PER  NODE 2 

MATERIAL  INDICATOR...  0 
GEOMETRY  INDICATOR...  0 
TWO-D  THEORY  INDIC...  0 

MAX.  NO  VERTICES 2 

MAX.  NO  WORDS/ELMT. . . 5 

NO  OF  DRC  UNITS A 

EXECUTION  OPTION A 

PRINTOUT  INDICATOR...  1 
RELABEL  INDICATOR....  0 
COORD  GEN  INDICATOR..  0 
DBG  GEN  INDICATOR....  0 
ELMT  GEN  INDICATOR...  0 
LOCAL  COOR  INDICATOR.  0 

Gt 

G2....... ............  . 

G3 

ACEL . 

NO.  OF  CONCENTRATED  LOADS. 

NO.  OF  PRESSURES 

NO.  OF  MATERIALS 

NO.  OF  THICKNESSES 

NO.  OF  TEMP.  INCREASES.... 

NO.  OF  Y-GRADI ENTS 

NO.  OF  ^.-GRADIENTS 

NO.  OF  CROSS  SECT.  AREAS.. 
NO.  OF  TORSIONAL  CONSTANTS 
NO.  OF  Y-MOM . OF  INERTIA.. 
NO.  OF  7,-MOM.  OF  INERTIA.. 


OF  LOCAL  COORD  SYSTEMS 


-MATERIAL  PROPERTI ES- 


MATL 

1 

2 


i 25000+06 
,50000+06 


G 

,00000 

,00000 


ALPHA 
. 00000 
. 00000 
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- KSH  POINT  DEFINITIONS 


NODE 

X 

Y 

7 ; 

1 

.0001) 

. 0000 

. 0000 

2 

.0000 

. 1 000+02 

. 0000 

3 

.0000 

.2000+02 

.0000 

A 

. 1000+02 

.1000+02 

.0000 

5 

. 1 000+02 

.2000+02 

.0000 

6 

. 1 500+02 

. 1 500+02 

. 0000 

7 

.2000+02 

. 1000+02 

. 0000 

8 

.2000+02 

.2000+02 

*0000 

9 

. 3000+02 

. 1000+02 

.0000 

10 

. 3000+02 

.2000+02 

.0000 

1 1 

. 3000+02 

.0000 

. 0000 

DEFLECTION  COWSTRA I MTS- 


DEPENDENT  INDEPENDENT 
NODE  DOF  NODE  DOF 


VALUE 


1 

I 

1 1 
11 


1 

2 

1 

2 


1 

1 

1 1 
1 1 


1 

2 

1 

2 


I 00000 
. 00000 
. 00000 
I 00000 


-ELEMENT  DESCRI PTI ONS- 


HMT  ELMT  MATL  NODE  NODE  ANOL  AREA  JMMX  JMMY  JMMX  PRFS  1FMP  .JSDY  .TSDT 
NO  TYPE  TYPE  l 2 TY1  E TYPE  TYPF  TYPE  TYPE  TYFE  TYPF  TYPE  TYPE 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


12  0 1 
2 3 0 1 

1 A 0 1 

2 A 0 1 

3 A 0 1 

3 5 0 2 

A 5 0 2 

A 7 0 2 

A 6 0 2 

2 6 0 2 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

c 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

n 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


5-15 


I 


1 1 

1 

2 

6 

8 

n 

1 

0 

0 

c 

0 

0 

n 

0 

12 

1 

2 

6 

7 

n 

1 

0 

0 

0 

0 

0 

0 

0 

13 

1 

P. 

5 

« 

0 

1 

0 

0 

0 

0 

0 

n 

0 

1 4 

1 

2 

7 

8 

0 

1 

0 

0 

n 

0 

n 

ri 

n 

1 5 

1 

2 

7 

1 1 

0 

1 

0 

0 

n 

0 

0 

n 

n 

16 

1 

2 

7 

9 

0 

3 

0 

0 

0 

0 

n 

0 

0 

17 

1 

2 

7 

1 0 

0 

3 

0 

0 

0 

0 

0 

0 

n 

18 

1 

2 

8 

10 

n 

3 

0 

0 

0 

0 

0 

n 

0 

13 

1 

2 

9 

1 1 

0 

3 

0 

0 

0 

n 

n 

n 

n 

20 

1 

2 

9 

10 

0 

3 

0 

0 

0 

0 

0 

n 

n 

-C  3NCENT RATED  MODAL  LOAT- 


NODE  DGF  VALUE 
5 2 - . 1 000+02 


EM;  r.F  input  for  linear  elastic  problem 


IS  THi'iE  ANOTHER  DATA  CASE  < YES  OH  JO 

:jo 


NORMAL  EXIT.  EXECUTION  TIME:  22^8  M ILLI  St:CONDS. 

DATA  10 NOR ED  - IN  CONTROL  MODE 


The  output  file  PR0BO1.  will  now  be  edited  to  change  the  area  information 
that  was  misinterpreted  by  the  program.  We  also  noted  during  the  summary  of 
input  data  that  the  connectivity  information  for  member  10  is  incorrect.  The 
element  is  presently  connected  between  nodes  2 and  6;  this  will  be  corrected 
to  be  between  nodes  5 and  6.  The  editor  session  to  correct  PR0BO1 . is  as 
follows. 
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-09: 50-c  nn): f 

kdit 
*o  4 

i .loo+on  a . aoo  o .oon 

* C /0  .000/3  .300/ 

i .loo+oo  a .auo  3 . 30c 

tX 

*P  00 


plane  thuss  prodlbm  i.o  manual  nh 


i i 20200025  4 1 

.l»O0 

l .800+ 06  .000  .non 

l .100  + 00  8 .'*00  .0  . 


0800003000  041000008) 


.0 


a . soo+oo  .ooo 


300 


1 

0000 

. oooo 

. oooo 

2 

oooo 

. 1 000+09 

. oooo 

3 

0000 

. 8000+0 

i *J 

. oooo 

/i 

1000+02 

. 1 000+02 

. oooo 

b 

1000+02 

. 2000+02 

. oooo 

6 

1500+02 

. 1 500+02 

.0000 

7 

8000+02 

. 1 000+02 

. oooo 

8 

2000+02 

. 2000+02 

.0000 

9 

3000+02 

. 1 000+02 

. oooo 

10 

3000+02 

.2000+02 

. oooo 

l 1 

3000+02 

.0000 

. oooo 

1 1 

1 1 

00 

12  12 

.00  111 

1 1 1 

- 1 

10  1 

100 

0 

1 8 

0 

0 0 

0 

0 

- 8 

101 

100 

0 

2 3 

0 

0 0 

0 

0 

-3 

10  1 

1 00 

0 

1 4 

0 

0 0 

0 

0 

- it 

101 

100 

0 

8 4 

0 

0 0 

0 

0 

-5 

101 

100 

0 

3 4 

0 

0 0 

0 

0 

- 6 

101 

800 

0 

3 5 

0 

0 0 

0 

0 

-7 

101 

200 

0 

4 5 

0 

0 0 

0 

0 

-H 

101 

200 

0 

4 7 

0 

0 0 

0 

0 

-9 

101 

800 

0 

4 6 

0 

0 0 

0 

0 

1 o 

101 

800 

0 

2 6 

0 

0 0 

0 

0 

1 1 

1 02 

100 

0 

6 8 

0 

0 0 

0 

0 

1 2 

102 

100 

0 

6 7 

0 

0 0 

0 

0 

1 3 

102 

100 

0 

5 « 

0 

0 0 

n 

0 

14 

102 

1 00 

0 

7 8 

0 

0 0 

0 

0 

1 5 

1 02 

1 00 

0 

7 1 1 

0 

0 0 

0 

0 

1 6 

102 

300 

0 

7 9 

0 

0 0 

f! 

0 

17 

102 

300 

0 

7 10 

0 

0 0 

n 

0 

18 

102 

300 

0 

H 10 

0 

0 0 

D 

0 

19 

1 02 

300 

0 

9 1 1 

0 

0 0 

0 

0 

20 

102 

300 

0 

9 10 

0 

0 0 

0 

0 

52 

“ • 

1000+02 

82  -.1000+02 

, 000 


00 


118  118 


OLND 

EPF  AT  LIME  38 

♦ L - 10 

-10  101  200  0 

♦ C /2  6/5  6/ 

0 


2 

5 


6 

6 


-10  101  200 

♦ E 

END  EDIT  38  LINES  OUTPUT 


0 0 0 0 0 


0 0 0 0 0 


.0 


. 00 
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5. 3. 1.2  Data  File  Node 

The  input  data  for  the  VISCEL/ELAS  IDP  is  created  by  using  the  text  editor. 
First  a file  must  be  created  using  the  system  commands 

(»CAT,UP  ELASPROBS 
0FREE  ELASPROBS. 

An  element  called  TRUSS  will  now  be  created  using  the  UoM  text  editor  as 
follows: 

@ED  TRUSS 
INPUT 

(The  input  data  has  not  been  shown) 


wmmmmm 


I 

I 

^ Note  that  the  element  TRUSS  was  created  In  the  workspace  and  copied  to  the 

file  ELASPROBS .TRUSS  by  the  command 

<»C0PY,  S TRUSS,  ELASPROBS. 

The  I DP  Is  now  executed  In  the  data  file  mode  as  follows,  where,  again  we 
assume  that  PR0BO1  contains  the  formatted  output  file. 


9XQT  symposium^-preprocessor.  v~  scel 

*****  TS007  *** 

I Tv  PUT  DATA  PREPROCESSOR  FOR  ELAS/VISCEL 

ENTER  NAME  OF  OUTPUT  DATA  FILE 
PROBOl. 


DO  /OU  WANT  TO  ENTER  DATA  FROM  A DATA  FILE  CYES  OR  MO) 
YF.S 


EM TER  NAME  OF  INPUT  DATA  FILE 
ELASPROBS. TRUSS 


-SUMMARY  OF  IT.  AS/ VI  SC  EL  INPUT  DATA- 


PLANAR  TRUSS  CASE 


-CONTROL  I NFCHMAT I ON- 


NUMBER  OF  NODES 11 

NUMBER  OF  ELt-lTS 20 

D.C.F  PER  NODE 2 

MATERIAL  INDICATOR...  0 
GEOMETRY  INDICATOR...  0 
TWO-D  THEORY  IMDIC...  0 

MAX.  NO  VERTICES 2 

MAX.  NO  WORDS/ELMT. . . 5 

NO  OF  DBC  UNITS 4 
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EXECUTION  GPTIGN 4 

PRINTOUT  INDICATOR...  2 
RELABEL  INDICATOR....  0 
COORD  GEN  INDICATOR. . 0 

DBC  GEM  INDICATOR....  0 
EL: IT  GEN  INDICATOR.  • • f) 
LOCAL  COOR  INDICATOR.  0 


G 1 .0000 

na .oooo 

G3 .0000 

AC  1 L .0000 

NO.  OF  CONCENTRATED  LOADS..  g 

NO.  OF  PRESSURES 0 

NO.  OF  MATERIALS P 

NO.  OF  THICKNESSES 0 

NO.  OF  TEI1P.  INCREASES 0 

NO*  OF  Y-G RA  i * I EN  TS. .•..«•*.  0 

NO.  OF  X-GRALT  ENTS 0 

NO.  OF  CROSS  SECT.  AREAS...  3 
NO.  OF  TORSIONAL  CONSTANTS.  0 
N* » • Oh  Y—  M OT  . OF  INERTIA...  0 
Ji  i . i jl*  ^ — M OT'l  . i « F I N E.AT  t A . . . 0 

NO.  OF  LOCAL  COORD  SYSTEMS.  C 


MATERIAL  PROPERTILS- 


MATL  E G 

1 • 850l'0+06  . 00°00 

g • 50000+  0 6 .00000 


-CROSS  c ECT  I ONAi . AHKAS- 


TYPL  VALI  E 
1 .10000+00 
g .aoooo+oo 

3 .30000+00 


-MESH  POINT  DEFINITIONS- 


NOLL 

X 

Y 

1 

. 0000 

.0000 

2 

. 0000 

. 1 000+02 

3 

.0000 

. 2000+02 

A 

. 1000+02 

. 1 000+02 

5 

. i ono+op 

.2000+02 

6 

. 1500+02 

. 1 500+02 

7 

.2000+02 

. 1 000+02 

8 

. 2000+02 

.2000+02 

9 

. 3000+02 

. 1000+02 

10 

•3000+08 

. 2000+02 

1 1 

. 3000+08 

.0000 

ALPHA 
. 00000 
. 00000 


. 0000 
.0000 
. 0000 
.0000 
.0000 
.0000 
. 0000 
. 0000 
.0000 
.0000 
. 0000 
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- I > EFL  EC  T I OM  C ON  STHAI  N T $- 


DEPENDENT 

INDEPENDENT 

NODE 

DOF 

NODE 

DDF 

VAL 

1 

1 

1 

1 

.00000 

1 

2 

1 

2 

. 00000 

1 1 

1 

1 1 

1 

. 00000 

1 1 

2 

1 1 

2 

. 00000 

EL EM EM T DESCRIPTIONS- 


MT 

ELM  T 

MATL  NODE 

NODE 

ANGL 

AREA 

•JMMX 

JMMY 

JMMX 

F*EC 

TEM  P 

jsr  ■> 

J c DZ 

0 

TYPE 

TYPE  1 

2 

TYFE 

TYPE 

TYPE 

TYFE 

TYPE 

TYPF 

TYPE 

TYFI- 

TY  &F 

1 

l 

1 

1 

2 

p 

1 

0 

0 

0 

r 

n 

n 

0 

2 

1 

1 

2 

3 

(i 

1 

0 

0 

0 

0 

0 

n 

0 

3 

1 

1 

1 

4 

0 

1 

0 

0 

0 

0 

0 

c 

0 

/i 

1 

1 

2 

4 

0 

1 

0 

0 

0 

0 

0 

r 

0 

5 

1 

1 

3 

4 

0 

1 

0 

0 

0 

0 

o 

0 

0 

6 

1 

1 

3 

5 

0 

2 

0 

0 

0 

0 

n 

0 

0 

7 

1 

1 

4 

5 

0 

2 

0 

0 

0 

c 

n 

n 

0 

8 

1 

1 

4 

7 

0 

2 

0 

n 

0 

0 

0 

0 

o 

9 

1 

1 

4 

6 

0 

2 

0 

0 

0 

0 

0 

a 

n 

10 

1 

1 

5 

6 

0 

2 

0 

0 

o 

0 

0 

0 

0 

1 1 

1 

2 

6 

8 

0 

1 

o 

0 

0 

0 

0 

0 

0 

12 

1 

2 

6 

7 

0 

1 

n 

0 

0 

0 

0 

n 

0 

13 

1 

2 

5 

8 

0 

1 

0 

0 

0 

0 

n 

0 

0 

14 

1 

2 

7 

8 

0 

1 

0 

0 

0 

0 

0 

0 

n 

15 

1 

2 

7 

1 1 

0 

1 

0 

0 

0 

0 

0 

n 

n 

16 

1 

2 

7 

9 

0 

3 

0 

0 

0 

0 

0 

0 

0 

17 

1 

2 

7 

10 

0 

3 

0 

0 

0 

0 

0 

n 

0 

18 

1 

2 

8 

10 

0 

3 

0 

0 

0 

0 

0 

c 

0 

19 

1 

2 

9 

1 1 

0 

3 

0 

0 

0 

0 

0 

0 

0 

20 

1 

2 

9 

1 0 

0 

3 

0 

0 

0 

0 

0 

0 

0 

C MM CENTRA  TED 

NO 

DAl»  L HAD  • 

NPRE  I >0F  VALUE 

a : -.iooo+02 


END  OF  INPUT  FOR  LINEAR  ELASTIC  PRGRLEM 


MORMAL  EXIT.  EXECUTION  TIME*  2771  MILLISECONDS. 
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5. 3. 1.3  Analysis  Problem  Execution 

Using  one  of  the  two  modes  in  the  IDP  we  have  created  a file  PR0BO1 . which 
contains  the  input  data  for  the  VISCEL/ELAS  analysis  program.  The  analysis 
program  can  be  executed  in  the  demand  mode,  but  since  the  output  may  be  more 
than  we  care  to  view  on  the  terminal  we  will  redirect  the  output  by  the  command 
^SUSPEND  and  resume  terminal  printing  by  PRESUME.  These  two  commands  are  trans- 
parent, and  serve  the  function  of  BRKPT'ing.  We  also  note  that  since  the  input 
is  in  file  PR0BO1  we  must  append  it  to  the  input  stream  by  an  @ADD-card. 


mpm^ si  ' -•i*psuGF.-;<A:-i . vi  scel 

Q ADI  FRGPOl. 

I.  IKS'.  ME 

LXAMINE,  PRINT*  DELETE,  GR  HOLD?  E 
EDIT 


The  results  will  be  scanned  to  determine  whetheA  they  agKee  with 
published  results^ the  intenative  session  is  as  follows. 


*L  NGDAL  D 
MODAL  DEFLECTIONS 
» P 10 

MODAL  DEFLECT I CMS 


GD” 

X-DISPLMT 

Y -DI  SPLMT 

-DI  SPLMT 

X -RGTATM . 

Y-RGTATM. 

7- RCTA 

1 

.0000 

.0000 

. 0000 

.0000 

. 0000 

.0000 

2 

“•636  6_  03 

-. t 566-02 

.0  100 

.0000 

.0000  ' 

.0000 

3 

-.  1954-02 

-.3133-02 

.0000 

.0000 

. 0000 

.0000 

A 

-.6366-03 

- .6247-02 

.0000 

. 0000 

. 0000 

.0000 

5 

-.  37  37-02 

-.7706-02 

.0000 

. 0000 

. 0000 

. 0000 

6 

-.  1721-02 

-.5926-02 

.0000 

. 0000 

.0000 

.0000 

* E 
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The  aeAuttA  agx.ee,  the  contenti  the  output  fate  axe  thu&  pointed 
at  RMENGR 


EXAMINE#  PRINT#  DELETE,  OR  HOLD?  _P_ 
WHERE?  RM ENG R 
“ SKAFI 


SENT  BY! 


■FT  S Ri-IFNGR 


5. 3. 1.4  Printed  Results 

The  output  file  for  the  example  problem  is  as  follows: 


W) 


tlam  -nu ee  prc-l-m  in  manual  nc  8 


?r"PLE«  -y~e:  equilibrium 


-C-AL  NUMEER  OF  NOTES  11 
-C-AL  NUME  rn  or  r LEME  NTS  20 
5Zrr>r~  Z ,'r  rRErEON  PER  NOC—  “* 
I T v F 0 VALUE  0 
I6'"M  VALUE  0 
ISTC  VALUE  0 
MO.Cr  V^RTTO'S  PER  ELEMENT  2 
V.'  C ° E S PER  ELEME  N-  "ESCPPTN.  E 
NO. O'*  0°“  TN?UT  UNITS  4 
NC  • Cr  C.LOAE  INPUT  UNITE  2 
NO  .Cc  "R^SSU^E  T YP1"  E 0 
NC.CF  MATERIAL  -VFEE  2 
NO.Cp  THICKNESS  TY°ES  0 
NC. CP  7E**P  .CHANCE  ’YPES  0 
NO. OF  Y-TMP.G9A9.  TYPES  0 
NC.CF  I--rEY°.CRAD  ,-YPES  0 
NO  •nr  A Rr  A T Y PE S 1 
NC • Cr  TORS-CN  CCNE.-YFES  0 
NO. OF  Y,-*'OMNT.OF  IN.TY0ES  0 
NC.CF  Z-MOKNT.OF  IN. TYPES  0 
NO. OF  ANCLE  FT  TYP"S  0 
’N*  V ?LUr  4 
IN”  V A L U " 1 
~EWLr  VALUE  0 
ICOP  VALUE  0 
-FUN  VALUE  0 
IM"S  VALUE  0 
IF-P  VALUE  0 
SCRATCH  TAPE  MUMPER  9 
nLOT  IN^UT  TAPE  NUM5ER  0 


ACrELEPATN*UNIT  v A EE  .OCCC 
ACPEL. OIR.COS.  .000  .000  .000 
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NC"«t  DEFLrCTCN': 


•'r  < - DIC^LMT 

Y 

-070PLMT 

Z- 0T5PLM  T 

X-ROTATN. 

y-->OT  AIN. 

7-ROT  AT 

T "(10'' 

r*  **rs  r* 

• L k_ 

.'’OCC 

.COCO 

• TCC 

..CTO 

2 - .3T6r-23 

- 

.1 ‘56-02 

• 0200 

.00  00 

»-*  o o o 

.G',T 

“•  _ i o c (i  _ n o 

• * _ — <_ 

- 

• 7 13 ’ -C  2 

• COCO 

.coco 

.^OCC 

..0',rC 

4 - .^Gr-?3 

- 

.5’’47-02 

.0000 

9 r*n  op 

.0"0" 

o o n r* 

• o 

r _ o-»7  -r  „nr» 

- 

.7700-02 

• 7000 

.coco 

.rocc 

..C',rC 

" -.1^21-?? 

- 

.5726-02 

.0000 

. C"OC 

.o^co 

. c n 

- .^999  -93 

- 

. 3072-02 

.0000 

.coco 

.'TCC 

. c-'-r 

r _ 7(?pr  — 

- 

.5431-02 

. 0000 

• 0 OC  0 

.rocc 

.C~rC 

n _ c ->  Qa.n  3 

- 

.2'll-03 

.0000 

.0^00 

# n r nn 

.C7'*n 

?T  -.3241-02 

- 

• S**2  1-D  3 

.0000 

.0000 

.rccc 

..C'TC 

ii  .2^0r 

^ o o n £ 

. 0000 

. 00  CO 

0 0^0 

• k..  ^ 

.-non 
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T'ro:  ant 

VCVEN7  "’rSULT  ants  at  t 

HO  NC or 

Cr  CN'  OTMCNETCML 

ELCMCNTC. 

uant:  ti-s 

A” 

E IN  LOCAL  COO'OIN  T 

0 3Y3TC 

M 

3 A'  TH'  ENOS  OF  ELEMFNT2 

L7  NCD  f r 

?;-x 

2-  Y 

f>  _ 

•y 

L. 

"-X  M - 

y v - 7 

^ P P 

y7 

.00  7 1. 

nop  * <">  n 0 •“  t 

-1.00° 

.coo  . 'pozt  .ccc 

• Cu  u 1 • D ° ^ 

1 11 

- 

. Z31 C +'l 

.crco 

.or  00 

.T3C  .C'CC 

.0000 

1 ? 1 

- 

.’OlE+Ol 

r pop 

.0070 

• oo"  .coco; 

.0027 

"JO  “CS 

y*r 

n r-  r -i 

• v-  1 • 

oco  .'cc:7 

-1. 200 

.00'  . 0 CC  2 7 .'O' 

.CCC  l.CCC 

0 7 1 

- 

. "312+01 

^ p 0 0 n 

# n n p 

.70"  .coco: 

.0007 

' 7 1 

- 

. 73 1 C +' 1 

.0000 

poop 

• l — L 

.'ore  .occc 

.CCC  c 

7"I'  "3 

XT 

■yr*  -+ 

• • 

7C  7 .0C0ET 

- .7C7 

.70,  , f 7'Z  T .OCT" 

.aoc  - .77' 

7 1 1 

- 

#orA»^4r*2 

P ^ P P 

• w c.  L.  - 

.”000 

.'or:  .o'cc 

.core 

7 A 1 

- 

. 7 0 ' 4 +0 1 

.cooo 

.0000 

.'CCC  .C'CC 

.core 

a"!"'  "S 

Y ▼ 

*3  nn  “» 

00"  .OCCOT 

# p no 

1 

.000  . CC'ZT  . CO' 

.002  1.0" 

«*  ° 1 

- 

#i  r.py 

.coco 

.0000 

,'rrc  .ccr.c 

.cocc 

«i  4 1 

- 

. 14  °7- 07 

.000" 

.0070 

.7  o"  .coca 

.0000 

"IF  -2r 

yT 

t r 7 

• • 

707  .'■'COOT 

.707 

. 7' 7 .7CCZT  .'C' 

.ccc  l.rcc 

- 7 1 

r r to  + ^ i 

."00 

.0070 

.70"  .coca 

^ pppn 

r 4 1 

•’rn?+rl 

.OCO' 

. '000 

.'ccc  .croc 

.cocc 

r-r>  -r-i  f*  f\  r 

XT 

1.007  . 

000  . 00  2ET 

.00” 

y 

X 

.000  • 'O'ZT  .COT1 

.aoc  1.7" 

' 7 1 

- 

• 731T  *01 

.occc 

. nocc 

.'err  .croc 

.oooc 

' 5 1 

- 

. "10  + 01 

O P ori 

• U ~ J -J 

.0000 

# p 0p  p # CC  DC 

An  rt  a 

• 0 L*  _ 

•"  1 

A.  w _» 

X” 

.007  1. 

O r-»  0 # n n Q r ^ 

-1.00' 

.000  ,'C'ZT  .00' 

.002  1.2'C 

7 4 1 

- 

,7230+'l 

.000' 

.occc 

.'CCC  .occc 

r-  r a r 

• w w < C 

- c 1 

- 

^ 7 2 o o + p 2 

• ccc  ^ 

. '000 

.'C'C  .occc 

.cc:c 

orx  — 

XT 

1.000  . 

-*no  nnr"T 

# 0 on 

1 

.007  ."'ZT  . 0C7 

^ p 

• Jju  X • . 

”41 

. r377+C0 

^ poor 

.0000 

.00"  .cooo: 

.coco 

”71 

. 3277  + 7C 

.000' 

.0000 

.oon'  .coca 

n +a  n ^4 
touju 

Of>  J-  ~ />  C 

X” 

7 7 . 00  0ET 

-.707 

. 707  . CC'ZT  .00' 

• aoc  1.0" 

o 4 1 

- 

.^s?r*ci 

.0000 

."CCC 

.'ccc  .core 

.C'CC 

' 7 1 

- 

-*  o op  4.  r*  2 

.0000 

.0000 

.000'  . 0020 

.027' 

2 **  *”  ▼ £ r <>  (* 

X” 

.7"  -. 

7 c7  .00C3T 

.70  7 

• 7'7  .CCCZT  .'CO 

.022  l.'CC 

r ^ 

- 

.zcor+ci 

.coco 

. oooc 

.'ere  .core 

.cocc 

in  si 

- 

. T270+C1 

.0000 

.0000 

.00"  .occc 

.0200 

uriF  scs 

XT 

. 7 T7  . 

7 07  ."COOT 

-.70  7 

. 7' 7 .CCCZT  .CC' 

.CCC  l.'CC 

l7  SI 

- 

“*  9 pp  ♦ 02 

0 0 0 n 

• - u U w. 

.0000 

.oo7'  .ooca 

.200' 

11  4 1 

- 

. ?370+"l 

.0007 

.0000 

.oc"  . o' on: 

.0020 

1"JF  70S 

y- 

7r^  _ 

• ’ • 

7 07  . TOO  7T 

.70  7 

. 7" 7 .CCCZT  .OC' 

.occ  i.rcc 

1”  si 

- 

# 7 0 02 

• COD? 

.0000 

.oc"  .ooca 

• no 

1”  7 1 

- 

.7E0t+01 

.coco 

.'COO 

.'ccc  .core 

.core 

2 **n  ~ r>r 

XT 

1.000  . 

00°  .00C0T 

.007 

1 

.00'  • "7ZT  .20' 

.002  ^.C" 

1’  3 1 

- 

.izi3.+n 

.coco 

.'CCC 

.'ore  .core 

.cocc 

1’  3 1 

- 

. 1213+01 

.0000 

.0000 

.'0"  .oroc 

. 0000 

: r4  r t r r £ c 

X7 

• ”00  1. 

CCC  .''COOT 

-1 .CCC 

.00'  .CCCZT 

.ccc  i-.rrr 

14  7 1 

- 

,7203+'l 

.occc 

.'OCO 

.'ere  .core 

.cocc 

14  3 1 

- 

,7293+71 

.0000 

.0000 

.0007  .C'CC 

.0200 

1"IF  OCO 

X7 

• 7 0 7 

707  ."COOT 

.70  7 

. 7 ' 7 .CCCZT  .'0' 

.OCO  l.CCC 

1-  7 1 

- 

• ® 0 °4  4-  PI 

.0000 

.0000 

.00"  .cocc 

.2000 

ir  :i  i 

- 

#°C'4+C1 

.coco 

.occo 

.'CCC  .O'CC 

.oooc 

l'OI'  rs 

XT 

1.000  . 

000  .0C00T 

.00' 

1 

.COO  . '0°  Z T .00' 

.000  1.0" 

i"  7 l 

- 

. II77- 06 

.0000 

.0000 

.'OCO  .0000 

.cooo 

lr  3 1 

- 

•1100-n6 

.coco 

.0000 

.'OOC  .cocc 

.oooc 

2 “'pi  r ^£0 

XT 

“*  r-  7 

• • 

707  • 'CCC? 

-.707 

• 7 ' 7 .CCCZT  .r.cc 

.occ  l.rcc 

1’  7 1 

.3373+01 

.0000 

.0000 

.7 on'  , .ooca 

.000' 

r to  i 

. 53T3+C1 

.0000 

.0000 

.oco'  .coca 

.0000 

1 ”' 17  "S 

X” 

1.000  . 

7oo  .ncocT 

.007 

2 

.000  .'O'ZT  .02' 

.aoc  l.rcc 

1°  3 1 

- 

.toic+'i 

.coco 

.occc 

.'ccc  .ccc.: 

.cocc 

1”  1 2 1 

- 

. 7013 ♦ "1 

^ n p p p 

.0000 

.OC"  .000" 

.coco 
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Table  5.1  - VISCEL/ELAS  Preprocessor 
Input  Data  Elements 

~ ■—  ' ^ 

Data  Item 

i. 

TITLE 

2. 

(IN.IT.IDEG) 

3.* 

(ITYPE,IGEM,ISTR,NVERT,NWDS,NDBC) 

4** 

(IMAT,NTIC,IAREA,IT0R,IYY,IZZ,IMFI) 

5. 

( I P , I PRS , ITEMP , I GRADY , I GRADZ ) 

6. 

LINK  - link  after  which  execution  is  terminated 

7. 

IPRINT  - print  parameter 

8. 

ILABEL  - relabelling  parameter 

9. 

. * 

ISHELL  - coordinate  selection  indicator  for  shells 

10. 

(G1,G2,G3) 

11. 

ACCELERATION 

12. 

Enter  IMAT  sets  according  to  value  of  ITYPE 

. EQ.  0 (E,G,a) 

. EQ.  1 (D]j » D12»  °22»  °24’  °44’  °55’  °56’  °66’ 

• EQ*  2 (°11*  °12’  °13*  °14’  °15’  °16’  D22’  °23’  °24’ 

0],  a2) 
°25’  °26 

°33'  °34’  °35*  °36’  D44’  °45’  °46*  °55’  °56’ 
°2*  °3} 

D66*  V 

13. 

IF(IPRS.EQ.O)  enter  no  data,  otherwise  enter  IPRS  pairs: 

(ISID,  PRESS) 

where  ISID  is  the  set  identification  number 

14. 

IF(NTIC.EQ.O)  enter  no  data;  otherwise  enter  NTIC  pairs 

( ISID, THICKNESS) 

5-28 

1 

J 

I 

I 

TABLE  5-1  (continued) 

^ 15.  IF(ISDT.EQ.O)  enter  no  data;  otherwise  enter  ISDT  pairs 

| (I SID,  TEMPERATURE) 

16.  IF(ISDY.EQ.O)  enter  no  data;  otherwise  enter  ISDY  pairs 

| ( ISID, Y-TEMP  GRADIENT) 

17.  IF  (ISDZ.EQ.O)  enter  no  data;  otherwise  enter  ISDZ  pairs 

(ISID,  Z-TEMP  GRADIENT) 

18.  I F( I ARE . EQ .0)  enter  no  data;  otherwise  enter  IARE  pairs 

(ISID, AREA) 

19.  IF(IMMX.EQ.O)  enter  no  data;  otherwise  enter  IMMX  pai'S 

(ISID,  TORSIONAL  CONSTANT) 

20.  IF( IMMY.EQ.O)  enter  no  data;  otherwise  enter  IMMY  pairs 

(ISID,  Y-MOMENT  OF  INERTIA) 

21.  IF;  “IZ.EQ.D)  enter  no  data;  otherwise  enter  IMMZ  pairs 

(IS10.Z-M0MENT  OF  INERTIA) 

22.  IF( IMFI .EQ.T)  enter  no  data;  otherwise  enter  IMFI  pairs 

(ISID, PHI) 

23.  ANS  (YES  OR  NO)  - IS  THIS  A PLANAR  PROBLEM? 

Enter  IN  data  sets 
. EQ.YES  (X,Y) 

. EQ.NO  (X.Y.Z) 

24.  If  (IBUN.EQ.O)  enter  no  data;  otherwise  enter  IBUN  sets  (NODE, I DOF,  NODE', 

IDOF' , VALUE) 

25.  ICLASS  - class  of  structure 

I 
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TABLE  5-1  (continued) 


25. A. 1 

25.B.1 
25. C. 1 
25. C. 2 

f 

I 

25.0. 1 
25. E. 1 

t 

25.F.1 
25. F. 2 

25. F. 3 

25. F. 4 

25  .G.l 
25.H.1 
25  .H.2 

| 

[ 

[ { 

* 


.EQ.l;  enter  NELM  - number  of  type  1 elements 
enter  NELM  data  sets 
(N1 ,N2 ,NMAT .IAREA, IPRESS, ITEMP) 

. EQ .2 ; same  as  26. A 

..EQ.3;  enter  NELM  - number  of  type  2 elements 
; enter  NELM  data  sets 
(N1 ,N2 ,NMAT , IAREA, I Z, IPRESS, ITEMP, I GRADY) 

. EQ .4;  enter  NELM  - number  of  type  4 elements 
; enter  NELM  data  sets 

(N1.N2.NMAT, IAREA, IY,IZ,IT0R,IANGLE, IPRESS, ITEMP, IGRADY.IGRADZ) 

.EQ.5;  enter  NELM  - number  of  type  3 elements 
; enter  NELM  data  sets 
(N1 ,N2,NMAT,IY,IT®R, IPRESS, IGRADZ) 

.EQ.6;  enter  NELM1  number  of  type  5 elements 
; enter  NELM1  data  sets 
( N1  ,N2 , N3 , NMAT , ITKN , I PRESS , ITEMP , I GRADZ ) 

IF(NELM1 . GE  IT)  enter  no  data;  otherwise  enter  NELM2  - number  of  type 
6 elements 

; enter  NELM2  data  sets 
(N1,N2,N3,N4,NMAT,  ITKN,  IPRESS,  ITEMP,  IGRADZ) 

. EQ. 7 ; same  as  27. F 

.EQ.8;  enter  NELMT1  - number  of  type  7 elements 
; enter  NELMT1  data  sets 
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TABLE  5-1  (continued) 


25.H.3 

25k  H.  4 

25k  1. 1 
25.1.2 

25. 1 . 3 

25.1.4 

25.J.1 

25.  J.  2 
25.J.3 

25. J. 4 

25.  K.  1 
25.  K.  2 

25.  K.  3 

25.  K.  4 


(ill  ,M2 ,N3, MMAT , ITKN, IPRESS, I GRADZ) 

I F( MELMT 1 . GE . IT ) enter  no  data;  otherwise  enter  NELMT2  - number  of 
type  3 elements 

; enter  MELMT 2 data  sets 
(Ml, M2 , M3 , M4 , MMAT , IT  KM .IPRESS.IGRADZ) 

. EQ.9;  enter  MELM1  - number  of  type  9 elements 
; enter  MELM1  data  sets 
(Ml ,M2 ,N3,M4,NMAT, IPRESS, ITEMP) 

; IF(MELMl.GE.IT)  enter  no  data;  otherwise  enter  NELM2  - number 
of  type  10  elements 

; Enter  NELM2  data  sets 
(Ml, M2, M3, N4,M5,M6,N7,M3, MMAT, IPRESS, MTEMP) 

.EQ.10;  enter  NELM1  - number  of  type  11  elements 
; enter  MELM1  data  sets 
(Ml ,N2, M3 , MMAT , ITKN , IPRE  SS , ITEMP , I GRADZ ) 

; IF(MELMl.GE.IT)  enter  no  data 
otherwise  enter  MELM2  - the  number  of  type  12  elements 
; enter  MELM2  data  sets 
(Ml, M2,. M3, M4.NMAT, ITKN, IPRESS, ITEMP, IGRADZ) 

.EQ.ll;  enter  NELM1  - number  of  type  13  elements 
; enter  .MELM1  data  sets 
(Nl, M2, N3.MMAT, ITKN, IPRESS, ITEMP) 

; IF(NELMl.GE.IT)  enter  no  data,  otherwise  enter  MELM2  - 
the  number  of  type  14  elements 
; enter  MELM2  data  sets 
(N1,N2,N3,N4,  MMAT, ITKN, IPRESS, ITEMP) 
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TABLE  5-1  (continued) 


25.L.1  .EQ.12;  enter  NELM1  - number  of  type  15  elements 

25. L. 2 ; enter  NELM1  data  sets 

(N1 ,N2,N3,NMAT,ITKN,IPRESS, ITEMP) 

25. L. 3 ; IF(NELM1 . GE . IT ) enter  no  data;  otherwise  enter  NELM2  - 

the  number  of  type  16  elements 
25.  L. 4 ; enter  NELM2  data  sets 

(N1 ,N2,N3,N4,NMAT,ITKN, I PRESS, ITEMP) 

25. M. 1 .EQ.13;  enter  NELM  - the  number  of  type  17  elements 

25. M. 2 ; enter  NELM  data  sets 

(N1,N2,NMAT,ITKN,IPRESS, ITEMP) 

25. N.1  ; enter  NELM  data  sets 

(N1 ,N2,NMAT,ITKN,IPRESS, ITEMP, IGRADZ) 

26.  IF(IP.EQ.O)  enter  no  data;  otherwise  enter  IP  data  sets 
(NODE, I DOF, VALUE) 

27.  ISUCA 

28.  IF(ISUCA.NE.O)  enter  3*.  4*,  followed  by  applicable  data  items  14-23  26-27. 
This  is  the  method  of  inputting  the  time  dependent  material  properties.  The 
termination  of  time  dependent  properties  is  indicated  by  setting  ISUCA  = 0. 

29.  ANS  - Response  to  question  "DO  YOU  WANT  A SUMMARY  OF  THE  INPUT  DATA  (YES 

OR  NO)" 

30.  ANS  - Response  to  question  "IS  THERE  ANOTHER  DATA  CASE  (YES  OR  NO)" 
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TABLE  5.3  Element  Properties 


6.  SATANS 


6.1  Introduction 

The  SATANS  program  solves  Sanders'  nonlinear  shell  equations  for  shells  of 
revolution  by  means  of  a finite-difference  approximation  in  both  spatial  and 

temporal  variables.  The  theoretical  development  for  the  program  is  described  by 

references  (9)  and  (10). 

The  governing  equations  are  reduced  to  coupled  sets  of  four  linear  second 
order  partial  differential  equations  in  terms  of  the  shell  meridional  coordinate 
and  the  time  by  expanding  the  surface  loads  and  the  dependent  variables  in  terms 
of  appropriate  fourier  series.  The  nonlinear  terms  are  treated  as  pseudo  loads. 

The  spatial  and  temporal  derivatives  are  approximated  by  finite  difference 
operator  having  orders  of  error  of  the  increment  squared.  A modification  of  the 
Houbolt  implicit  backward  difference  scheme  for  time  derivatives  is  used  to  insure 

stability  of  the  numerical  determination  of  the  time  response. 

The  applied  loads  may  consist  of  any  combination  of  pressure  loads,  thermal 
loads,  and  boundary  loads  and  initial  conditions  that  are  symmetrical  about  a 
datum  plane  in  the  shell. 

6.2  SATANS  Input  Data  Preprocessor 

The  SATANS  IDP  operates  in  the  demand  mode  and  allows  the  user  to  supply 
unformatted  input  data  in  a program-controlled  logical  sequence.  The  IDP  performs 
some  elementary  checks  and  writes  a formatted  input  file  for  the  SATANS  analysis 
program. 

The  IDP  operates  In  either  of  two  input  modes:  interactive  and  data  file 

modes.  In  the  interactive  mode  the  IDP  prompts  the  user  for  data  elements  which 
are  then  entered  from  the  terminal.  In  the  data  file  mode  the  user  pre-prepares 
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the  responses  to  the  IDP  input  requests  and  saves  them  on  a cataloged  file.  The 
IDP  will  then  by-pass  all  prompting;  and,  will  read  the  unformatted  input  from  the 
data  file. 

In  either  mode  the  IDP  will  print  a summary  of  the  SATANS  input.  This  summary 
is  actually  obtained  by  reading  the  formatted  SATANS  input  file  thus  providing  a 
check  on  both  the  correctness  of  the  data  items  and  the  file  structure. 

The  input  data  items  required  by  either  mode  are  the  same.  These  data  items 
are  presented  by  Table  6.1. 

6.2.1  Interactive  Mode 

The  IDP  is  executed  in  the  interactive  mode  by  responding  "NO"  to  the  query, 
"Do  you  want  to  enter  data  from  a data  file?"  The  IDP  will  then  lead  the  user 
through  the  creation  of  the  SATANS  input  file  by  prompting  for  data  elements  which 
the  user  is  to  enter  from  the  terminal  device. 

The  IDP  will  set  certain  parameters  by  default  according  the  users  specifica- 
tion of  a static  or  dynamic  analysis.  Thus,  for  example  the  parameters  that  are 
used  only  for  a dynamic  analysis  are  not  requested  if  ISORD  = 0. 

6.2.2  The  Data  File  Mode 

The  data  file  mode  allows  the  user  to  pre-create  the  responses  to  the  IDP 
requests  for  data  according  to  Table  6.1  and  thus  to  completely  by-pass  the 
prompting  mode.  In  either  the  data  file  or  interactive  modes,  the  user  may 
request  that  a summary  of  the  input  data  be  printed  at  the  terminal. 

The  new  user  may  find  it  convenient  to  use  the  IDP  in  the  full  prompting 
interactive  mode  as  a teaching  aid.  The  IDP  will  fully  describe  each  of  the 
data  elements,  and  will  provide  the  logical  branches  for  input  data. 
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the  associated  geometric  parameters  by  input  items  11  on  Table  6.1. 

The  explanation  of  the  input  parameters  requested  by  the  IDP  for  the  various 
shape  options  is  given  by  figures  6.1,  6.2  and  6.3  for  the  sphere,  cone  and 
cylinder,  respectively.  The  user  is  referred  to  the  users  manual,  reference  10 
for  a description  of  the  shell  parameters  for  the  case  of  a general  shell  of 
revolution. 

6.4  User  Prepared  Subroutines 

The  user  must  provide  subroutines  which  define  the  surface  pressure,  the 
thermal  loads,  the  stiffness  and  the  initial  condition.  The  function  of  these 
subroutines  is  described  below. 

6.4.1  Stiffness 


The  stiffness  parameters  are  calculated  by  the  subroutine: 


r 


The  nondimensional  stiffness  quantities  b,  db/d;,  d,  and  dd/df 
are  defined  in  BDB  for  each  meridional  station.  The  correspondence 
between  the  stiffness  quantities  at  the  Kth  station  and  the  F0RTRAN 
variables  is  as  follows : 

B - <b>K  = (BV<W 

= <§Vett) 

o o 

D = (d)K  = (D)K/(Eoh3) 

WK  'ds  K E h3; 

o o 


where 

B = JjEd£/(l-v2) 

D = Jc^dC/U-v2) 

6.4.2  Surface  Pressure 

The  surface  pressures  for  the  MNMAX  fourier  coefficients  are  to  be  defined 
by  the  subroutine: 

PLOAD  (K) 

The  nondimens ional  Fourier  coefficients  of  the  meridional,, 
circumferential  and  normal  components  of  the  pressure  load  p 'n' 

(n)  (n)  . s’ 

Pq  j and  p respectively , are  defined  in  PL0AD  for  each  meridional 

station  as  a function  of  the  Fourier  index. 

K(M)  . (P<”>)K-  (,s(n>)K  (a/ohj 


o o' 


pro.)  . (P<”>)K  . (,6'">)K  J 


o o' 


M = 1,  2,  ...  MNMAX 


ra(M)  = (P(n)L  = (q(n))„  (a/a  h ) 

iv  IV  O O 


Note  that  these  are  stored  as  functions  of  M only. 
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6.4.3.  Thermal  Load 


The  thermal  load  and  moment  coefficients  are  to  be  defined  by  the  subroutines 
TLOAD(K) 

The  nondimensional  Fourier  coefficients  of  the  thermal  loads 

trp^,  mT*'n\  ^ and  (m^/n  ) are  defined  in  TL0AD(K) 

for  each  meridional  station  as  a function  of  the  Fourier  index.  The 
F0F.TBAN  variables  are  defined  as  follows : 

TT(M1  = (tT(n))K=  (eT(n))K/(aoho) 

DT(M)  = (£  (t^)),  = (^)K(a/aoho) 

EMT(M)  = (mT(n))K  = (HT(n))K  (a/aoho3) 

Jn) 

DMT(M)  = (^(mT(n)))K  = (-Ir-)  (a2/coho3)- 

6.4.4  Initial  Conditions 

For  the  case  of  a dynamic  response  the  initial  conditions  are  defined  by  the 
subroutine: 


INITL 


Z(I,L)  = (z(n))K 


U(n)(Eo/aao) 

V(n)(Eo/aao) 

W(n)(Eo/aao) 

Msn)(a/ooho3)_ 


dz(n) 

ZD0T(I,L)  = ) 


- t — 
o dT 


U(n)(Eo/aao) 

v(n)(Eo/aao) 

w(n)(Vaao) 

^n)(Va0h03) 


I I = 1,2, 3, L 

L = 1,2,  ... 

(KMAX+2)*(MNMAX) 


The  index  L runs  from  1 to  KMAX+2  for  NN(l),  and  from  l+KMAX+2  to 
a(KMAX+2)  for  NN(2),  etc.  The  first  element  for  each  value  of  n 
corresponds  to  the  initial  ficticious  station,  the  next  element 
corresponds  to  the  first  station  on  the  shell,  etc. 
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6.4.5  COMMON  Data  Blocks 

The  parameters  required  to  define  the  variables  and  the 

variables  themselves  must  be  made  available  to  the  various  subroutines.  Generally 
speaking  these  data  elements  are  transmitted  by  NAMED  COMMON  statements.  The 
data  elements  and  the  associated  NAMED  COMMON  are  presented  by  Table  6.2. 

6.4.6  User  Definition  of  Subroutines 

Due  to  the  unusual  nature  of  this  program  the  user  must  modify  existing 
subroutines  prior  to  execution  of  the  SATANS  analysis  program.  The  user  thus 
proceeds  in  a radically  different  manner  to  execute  SATANS  than  is  used  to 
execute  the  other  programs.  Rather  than  using  the  normal  @XQT  statement  first 
proceed  by  making  the  following  system  command 

mO  SYMPOSIUM*PROGRAM.  SATANS 

This  statement  will  copy  several  elements  into  the  workspace  including  the 
following  default  symbolic  elements  for  user  prepared  subroutines: 


■ 


r 


i 

i 


7L  C 

c'  •••"*:  :/:cf  *4.r  j .r? (**.'>?:). :'iz: 

i /-.  1-4/ : t ( 4 . ?.zn r'\.r.r  ti '» »?2n)  *; pt * o~.t*  alca~ 

1/  • L /f.M  l"  ) . wf.TM' 
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A't  * J 
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" ” m y ■*. 

:/'l:  la :t  »ch’  k.  :r  c.-ibl'/4'*  ( nvminit 

?/~L7 ‘♦•r  *»r  > * s?  7sr»  al~  * r 
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7?  WC'J/3L3/P  ( 7 77  )»  ' VM  772)  . CV1T  ( ?~2> 

:/:-l'/l  : 7: r .it: 
i /•’l  * o l?  ?*?j 
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r 4 


A 7 { 1 ) 

jL?',C’JTI%C  7L0A7(X> 
r"tL  NU 
C?  Y*19N 

i/:r  Li/KNKAY/rsL'VNN  1 1 r ) 
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A4  9«» 

! 4 r 

A4  iZ 
A 42  7 
A423 

A4?° 
AE!  4 


A 4 7 1 
A4  7"* 
A473 
A4  74 
A 47  c 
A47S 
A 47  7 
A473 
A47*» 
A48? 
A 4 E 1 


A4  8T 
A 46  4 
A4er 
A 46  6 
A4  37 
A 46  ? 
A489 
A 427 
A491 
A 42  ? 
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The  user  may  now  modify  any  one  of  the  symbolic  elements  by  using  the 


1 


U of  M text  editor  (as  we  will  demonstrate  In  a later  section). 

6.4.7  SATANS  Analysis  Program  Execution 

Once  the  user  has  completely  defined  the  desired  subroutines  the  SATANS 
analysis  program  is  executed  by  the  command: 

PADD  XQTSATANS 

followed  by  the  addition  of  the  input  file  with  the  command 
PADD  FILENAME. 

As  in  the  other  analysis  systems,  the  user  may  wish  to  PSUSPEND  and  PRESUME. 

Let  us  suppose  that  the  IOP  had  created  the  input  file  SAT0UT  then  the  program 
would  be  executed  by  the  commands 

PSUSPEND 
PADD  XQTSATANS 
PADD  SAT0UT. 

PRESUME 

If  you  are  curious  about  what  the  element  XQTSATANS  Is  doing  you  may  list 
It  using  the  text  editor. 

6.5  Example  Problem 

Consider  the  problem  of  determining  the  linear  static  response  of  a cylindrical 
shell  subjected  to  the  edge  moments  as  shown. 
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(n)  = 


E = 30  x 106  psi 


v = . 3 

t = 1 inch 


56i* 


(n  = 0,  2,  5,  20) 


All  surface  pressures,  thermal  loads  and  initial  conditions  are  zero. 
6.5.1  Creation  of  SATANS  Input 

We  will  use  dimensional  input,  and  will  request  dimensional  output.  The 
interactive  execution  of  the  I0P  is  as  follows;  where  the  data  file  SATIN 
is  to  be  created. 
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•XQT  lYNPOSIUMPPRIPROCESSOR.  SATAN S 


•••••  PREPROCESSOR  FOR  SATANS 
00  YOU  WANT  MINIMUM  PROMPTING  (YES  OR  NG> 

NO 


INTER  NAME  OF  OUTPUT  DATA  FILE 

amour 


DC  YOU  WANT  TO  INTER  DATA  FROM  DATA  FILE  (YES  OR  NC> 

JS2 


INTER  PROBLEM  DESCRIPTION  (.LE«78  COLS) 
CYLINDRICAL  TEST  CASE 


•CONTROL  INFORMATION* 

INTER  THE  DATA  SET  ( ISORD#KMAX»MNHAX#NAXN> 

130RD- ANALYSIS  TYPE  INDICATOR 
•EQ*0*  STATIC  ANALYSIS 
• «E*U  DYNAMIC  ANALYSIS 

KM  AX  *N0  OF  EQUALLY  SPACED  MERIOIAL  STATIONS 
MNMAXRNO  OP  FOURIER  TERNS  USED  TO  DESCRIBE  B.C.  AND  LOADS 
MAXM  «N0  OF  FOURIER  TERMS  IN  SOLUTICM 
o»ae»A#4 


ENTER  A NODE  NUMBERS 


-STATIC  ANALYSIS- 

INTER  THE. DATA  SET  (LSMAXtLCHMAX#  I TRNAX#  DILOAD* EPS) 
LSNAXMNO  OF  LOAD  INTENSITIES 
LCHMAXrNC  OF  LOAD  INCRSMMT  REDUCTIONS 
1TMAXRMAX  no  of  itemticms 
dbaasrload  incrmbit 
EPS  ■ CONVERGENCE  CRITERION 


•PHYSICAL  PROPERTIES- 

ARE  PROPERTIES  (I*E.  LOADS  GEOM  ETC)  IN  DtNWSIONAL  PORN 
(YES>  OR  NO) 
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r 


I 


I 


I 


MTER  POISSONS  RATIO 


ENTER  THE. ORTA  SET  ( SI 00# BLAST* TXN* CHAR) 
S1G0  •RET*  STRESS 
■.AST* REP*  MODULUS  OF  SLASH  CITY 
TUN  •REF.  THICKNESS 
CHAR  •CHARACTERISTIC  SHILL  DIMM  SIGN 
I .»  1 »»  I »*  50* 


IS  THE  OUTPUT  TO  BE  IN  N CM  DIM  D)  SI  ONAL  FORM  (YES  OR  NO) 
YES 


- BOUNDARY  CONDI  TIMS' 

IS  SHELL.  CLOSED  AT  INITIAL  STATHM  (YES  OR  NO) 

JU- 
DO YOU  WISH  TO  MTER  ONLY  THE  DIAGONAL  ELMENTS  OF  CM  EGA  AND  LAM  DA 
MATRICES  (YES  OR  NO 
JUL 

MTER  THE  ELMENTS  OF  THE  MEGA  MATRIX  (BY  RONS)  AT  THE 
FIRST  STATION 
RON  1 

I «»0»*0«»0* 

ROW  2 

0.* 1**Q**0. 

ROW  S 

0**0** l**0. 

ROW  A 

0»*0*»Q.*0* 

MTER  THE  SL MINTS  OF  THE  LAM DA  MATRIX  (NY  ROWS)  AT  THE 
FIRST  STATION 
ROW  I 

A. *0**0. *0* 

ROW  fi 

o.*o**o*»o. 

ROW  3 

0.*  0**0. *0. 

ROW  A 

MTER  B.MMTS  OF  THE  ■.•VECTOR  AT  THE  FIRST  STATION 
0.*0.*0**.M 

IS  mmx  CLOSED  AT  FtNAL  STAT1CN  (YES  OR  NO) 

JUL 

BO  YOU  WISH  TO  MTER  ONLY  THE  DIAGONAL  N.MMTS  OF  OMEGA  AMD  LAM  DA 
MATRICES  (YES  OR  NO) 

IB 

MTER  DIAGONAL  B.MMTS  OF  HE  ONEGA  MATRIX  AT  THE  LAST  STATION 

°MTER  DIAGONAL  ELM  MTS  OF  THE  LAM  DA  MATRIX  AT  THE  LAST  STATICS) 
l/.0.*0*T 

f OtA-L* 

MTER  B.MMTS  OF  THE  ■.•VECTOR  AT  THE  LAST  STATION 
G.»t*»0.#.0S 
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•OUTPUT  CONTROL* 

HTU  THE  BATA  SET  (MINMAX* IFRSO* IPRIMT*  INODE) 
MOMAfiMMO  OP  ANGULAR  POSITIONS  < LE.6) 
inm  A PRINT  PREGUBfCY  ALONG  A NCR  I DICM 
I PRINT*  SOLUTION  PRINT  FABQUfiNCY 
INODE  -NODS  DATA  INDICATOR 

• EQ.O  SUMMED  OUTPUT 

• EG*  1 NODAL  OUTPUT 

ftlWUl. 


-GEOMETRIC  PARAMETERS- 
BITER  GEOMETRY  CODE  NUMBER 
•EG*  I SPHERE 

• EQ.fi  CONE 

• EG*  3 CYLINDER 
•EG. A GENERAL 


YOU  HAVE  SPECIFIED  CYLINDRICAL  GEOMETRY 
INTER  THE  DATA  SET  ( RADI  US#LBt6TH> 


BO  YOU  WANT  A SUNNARY  OF  INPUT  DATA  (YES  OR  NO> 
YES. 


SUNNARY  OF  SAT  MIS  INPUT 

CYLINDRICAL  TEST  CASE 


•CONTROL  INFORMATION- 

SQRD  A 0 
KHAX  A fifi 

mmaxa  a 

MUM  A A 
LSMAN*.  I 
LCM4AXR  I 

ITRNANA  I 

IC  A « 

TEEO  A •OGGOfi 

DILOAO  .10000*01 
MS  A • I60fi0*01 

I TAPE  • 0 
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•MY Si  CAL  PARAMETERS' 
HU  n . »»N«M 

si oo  a . toeoo-si 

ILAST  n •IOO0O+O1 

THICK  ■ .IMOO^Ol 

CHAR  • . SOOOO+Oi 


•OUTPUT  CONTROL  PARAMETERS' 
IMODE  ■ 1 

NDIMEM-  I 

NTHMAX-  0 

I FREQ  • 1 

I PRINT-  1 


•NODE  NUMBERS' 

THE  FOLLOWING  MODES  ARE  INCLUDED  IN  THE  ANALYSIS 
0 8 5 80 


-BOUNDARY  CONDITIONS' 
FIRST  STATION 


•ON EGA  MATRIX' 


•10900*01  *00899 

•00000 

•00000 

• 00000  • t 0099*01 

•00000 

•oooeo 

*00000  *00000 

•10000+01 

•09000 

•00000  *00999 

•09000 

.00000 

-LAM DA  MATRIX- 

•99999  .00009 

•90090 

•00000 

•99900  .00000 

•00000 

•00009 

•09000  .00000 

•00000 

•00099 

•09000  .00000 

•ootto 

•19909+01 

-EL  VECTOR- 

•00000  *00000 

•00000 

•90900-01 

LAST  STATION 
-ON E9A  MATRIX- 

•00000  .00000 

•00000 

•00090 

•00000  •! 0000+01 

•00000 

•99009 

•00090  .00009 

•10000+01 

•99999 

•09999  .90999 

•00000 

•99999 

•LAM DA  MATRIX- 

•19099+01  .90999 

•99999 

•99999 

•99999  .90999 

•99999 

•00999 

• 

3 

o 

• 

• 

e 

| 

• 

•90000 

•00999 

•00090  .00000 

•00000 

•19999+91 

-EL  VECTOR- 

• 

i 

• 

• 

• 

i 

•00000 

I 

• 

01 

• 
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•tCHITH!C  PARAMKTKRS- 


CYLlMDftlCAL  QE£MCTRY 

RADIUM  •5HI«*«a 
LDM1M* 


IS  1MCRC  AM  OTHER  CASK  (YES  OR  MO) 

HQ 


6.5.2  Listing  of  SATANS  Input  File 


The  SATANS  input  file  created  by  the  I DP  is  as  follows: 


7 

CD  fiOD  OS/06-lO«ie-(0O>tF 
EDIT 

♦y  so 

CYLINDRICAL.  TEST  CASE 


1 0 1 

Q 

10  11 

0 0 80 

4 4 

30000+00  • 

10000+01  *10000+01 

•10000+01 

•50000+00  • 

I 0000+01  • 

A 

10000+01 

V 

0 8 S 

00 

•10000000+01 

•00000000 

•00000000 

•00000000 

•00000000 

•10000000+01 

•00000000 

•00000000 

•00000000 

•00000000 

•10000000+01 

•00000000 

•00000000 

.00000000 

• 00000000 

•00000000 

•00000000 

•00000000 

•00000000 

•00000000 

i 

s 

• 

• 00000000 

.00000000 

•00000000 

•00000000 

•00000000 

•00000000 

•00000000 

•00000000 

• 10000000+01 

•00000000 

•00000000 

• 00000000 

•80000000-01 

•00000000 

•O0OO0000 

•00000000 

• 00000000 

•Wiffffl 

•10000000*01 

•00000000 

.00000000 

•00000000 

•00000000 

•10000000+01 

•00000000 

•00000000 

•00000000 

•00000000 

•00000000 

•10000000+01 

•00000000 

•00000000 

• 00000000 

•00000000 

•00000000 

•00000000 

• 

i 

o 

9 

•00000000 

•00000000 

.00000000 

• 00000000 

•00000000 

•10000000+01 

•00000000 

•00000000 

•00000000 

•80000000-01 

3 

•50000+00  • 50000+02 

tor  AT  LINE  80 
1L 

END  EDIT  M LINES  OUTPUT 
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6.5.3  Modification  of  User  Prepared  Subroutines 

The  default  symbolic  elements  for  the  subroutines  are  brought  into  the 
workspace  by  the  command 

@ADD  S YMPOSI UM*PROGRAM . SATANS . 

Considering  the  listings  of  the  symbolic  elements  presented  in  section 
6.4.6  we  see  that  only  the  stiffness  subroutine  BOB  must  be  changed. 

Since  the  input  is  dimensional  we  must  define  dimensional  quantities  for 
B & D which  are  constants  given  by 


D = 


Et' 


12(1  - v2) 


The  subroutine  BOB  is  thus  modified  as  follows  using  the  text  editor. 


ip  m 

CD  BSD  05/06-101 1 5-< CO)lF 
EDIT 


jt-fg 


SUBROUTINE  BDBC  K# B#  DB#  D»  DD1 
REM.  HU 

common 

l/M.3t/TKN*H#AST«CHAIUSI60/BLt5 

8 /W*  UK  19>#  V}<  !0>*Vtt  |0>»  V0<  1 0>#  UOC  I •)# Wtt  1 0>»U3<  10>»  V3<  I0>*  V3C  10) 
ypir/DiL 
RKTUm 

mo 

COT  AT  LIMC  8 

•u  a 

RETUIW 

•IL 

3/BL1T/DBL 

♦TAB  • 1 


6-15 


ILLEGAL  CHARACTER  SEQUENCE 
INPUT 

*D#B/1E« 

tPBfJLi — 

»C-PP»Qa 


EDIT 

liLJL 

tflJL 


■& 


D-B/12. 
DB*»0« 

EDIT  IP  * twrc 


6.5.4  Analysis  Program  Execution 

The  analysis  program  is  executed  using  the  data  file  SATANSIN.  as  input  as 
follows: 


•SUSPEND 
tAPD  XQTSATAHS 
•ADD  SATCWT. 

•RESUME 

examine*  print,  delete#  or  hold?  jl 

EDIT 
♦L  MS* 

•MAD  MAP#  SATAN  S 

•L_ 


•EOT 

SATAMSi 

•L  MS 

SPA 

MS  M TH  M ST 

Q S 

M S 

M T 

N ST 

•u 

STA 

MS  M TH  M ST 

A S 

M S 

M T 

N ST 

1 

• IM*M  #335“  D4  #tM 

-.167-02 

• IAA+91 

•30A*00 

• 000 

• 

-•3A4+03  *11144  •••# 

-•TOO- 04 

• S3S  «AI 

• TOT*  *8 

•AAA 

3 

••311*13  • |78*  §4  •Mt 

•• ITS* AS 

• 556“  03 

• I6T* 03 

• 000 

4 

••3lt*A3  *MM9  *111 

*•  4M*IT 

• 131*04 

• 39  3-OS 

• 000 

EXAMINE#  PRINT#  DELETE#  OR  KO.DT  JL 

VHERET  MflNR 

SUIT  BY!  SHAPES*  » EM— R 
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6.5.5  Analysis  Program  Output 
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Table  6.1  - SATANS  Preprocessor  Input  Data  Elements 


1 . 72  Column  Title 

2.  (ISORD,  KMAX,  MNMAX,  MAXM) 

3.  (NN  (I),  I ■ 1,  MAXM) 

4.  A If  (ISORD. EQ.O)  enter  following  data  set;  otherwise  go  to  4.B 

(LSMAX,  LCHMAX,  ITRMAX,  DELOAD,  EPS) 

4. B  If  (ISORD. GT.O)  enter  following  data  (LSMAX,  ITRMAX,  IC,  TEEO,  DELOAD, 

EPS,  ITAPE) 

5.  A Answer  (Yes  or  No)  to  query  "Are  properties  in  dimensional  form?" 

5.B  Enter  Poisson's  Ratio 

5. C  If  answer  Is  "no"  enter  the  following  data  set 

(SIGO,  ELAST,  TKN,  CHAR) 

6.  Answer  (Yes  or  No)  to  query:  "Is  the  output  to  be  in  nondimensional  form?" 

7.  A Answer  (Yes  or  No)  to  query:  "Is  the  shell  closed  at  the  first  station?" 

7.B  Answer  to  query  "Do  wish  to  enter  only  the  diagonal  elements  of  the  OMEGA 
and  LAMDA  matrices?"  (yes  or  no) 

7.CY  If  answer  is  "yes"  enter 
7.C.1  Diagonal  Elements  of  OMEGA 
7.C.2  Diagonal  Elements  of  LAMDA 
7.CN  If  answer  is  "no"  enter 
7.C.1  All  elements  of  OMEGA  (by  row) 

7.C.2  All  elements  of  LAMDA  (by  row) 

7. C.3  For  either  'yes'  or  'no'  enter  elements  of  EL-VECTOR 

8. B  Same  as  7.B 
8.C  Same  as  7.C 


L 
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Table  6.1  (Concluded) 


9.  (NTHMAX,  I FREQ,  IPRINT,  IM0DE) 

10.  If  ( NTHMAX. GT.Q)  enter  data  set  (ANGLE(I),  1=1,  NTHMAX) 

11.  Enter  GEOMETRY  CODE  (IGEOM) 

. EQ .1  Sphere 
. EQ .2  Cone 
. EQ.3  Cylinder 
.EQ.4  General  Shell 

11.1  If  (IGEOM. EQ.l) 

(PHIMIN,  PHIMAX,  RADIUS) 

11.2  If  (IGEOM. EQ. 2) 

(RMIN,  RMAX,  ANGLE) 

11.3  If  (IGEOM. EQ.3) 

(RADIUS,  LENGTH) 

11.4  If  (IGEOM. EQ.4) 

11. 4A  DEL  ^ 

11. 4B  KMAX  SETS  (R,  OMT,  OMXI,  DEOMX,  GAM) 

12.  Answer  to  query  "Do  you  want  a summary  of  the  Input  data?"  (Yes  or  no) 

13.  Answer  to  query  "Is  there  another  case?"  (Yes  or  no) 


6-27 


Table  6.2  COMMON  Data  Elements 


COMMON 

Data  Elements 

/BL  32/ 

TKN,  ELAST,  CHAR,  SIGO 

/BL 15/ 

NU 

/BL17) 

DEL 

/ I BL2/ 

NN(10) 

/ IBL1 ) 

MNMAX 

/BL3/ 

PR(10) , PX(10),  PT(10) 

/ IBL4/ 

KMAX,  KL 

/BL5/ 

TT(10),  EMT(IO) , DT(10),  DMT  (10) 

/BL101/ 

Z0(4 ,220) , Z2(4,220),  Z3(4,220) 

/BL1 00/ 

SORD,  TEEO 

/ IBL9/ 

MAXM 
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7.0  SAASIII 

7.1  Introduction 

SAASIII,  is  a general  purpose  finite  element  program  for  the  stress 
analysis  of  axisymmetric  and  plane  solids  with  different  orthotropic  tempera- 
turesdependent  material  properties  in  tension  and  compression.  The  program 
theoretical  and  user  information  are  presented  by  reference  (11). 

This  report  describes  the  execution  of 

(1)  An  interactive  input  data  preprocessor  which  produces  the  formatted 
input  data  file  for  the  SAASIII  analysis  program. 

(2)  The  execution  of  the  analysis  program  on  the  University  of 
Maryland  UN I VAC  1108. 

7.2  SAASIII 

The  SAAS  IDP  operates  in  the  Demand  mode  and  allows  the  user  to  input 

unformatted  data  in  a program-controlled  logical  sequence.  The  IDP  performs 
some  elementary  error  checks  and  writes  a file  which  contains  the  formatted 
input  modes,  interactive  or  data  file;  and,  further  allows  the  user  to  control 
the  amount  of  computer  generated  prompting  in  the  interactive  mode. 

In  the  interactive  mode  the  IDP  issues  prompting  information  which  is 
displayed  by  the  terminal  device.  The  user  then  inputs  the  requested  information 
directly  by  means  of  the  keyboard  in  unformatted  form;  i.e.,  the  data  elements 
are  delimited  by  commas.  This  interaction  continues  until  all  of  the  input 
requirements  have  been  satisfied.  In  the  data  file  mode  the  user  prepares  the 
unformatted  input  to  the  IDP;  and,  then  informs  the  IDP  the  name  of  the  data 
file.  The  IDP  then  reads  the  input  data  directly  from  the  data  file  in  exactly 
the  same  sequence  as  the  interactive  mode;  but,  completely  bypasses  the  prompting 
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In  either  mode  the  sequence  required  for  the  IDP  input  is  shown  by 
Table  7.1.  For  the  most  part  the  input  data  elements  have  the  same  names  and 
meaning  as  those  variables  defined  in  the  program  users  manual.  However  it  is 
recommended  that  the  interactive  mode  be  executed  initially  with  maximum  promp- 
ting in  order  to  completely  define  each  data  element. 

The  user  will  note  that  the  IDP  does  not  request  all  of  the  data  items 
described  by  the  user's  manual.  Specifically  the  IDP  does  not  request  any 
parameters  which  are  associated  with  plotting. 

7.2.1  Interactive  Mode 

In  the  interactive  mode  the  IFP  sets  up  a question-answer  dialogue 
with  the  user.  The  IDP  defines  a data  set  to  be  input  and  then  accepts  the  input 
in  unformatted  form  directly  from  the  terminal  keyboard.  The  SAASIII  IDP  has 
twe  levels  of  prompting  which  are  controlled  by  user  request.  The  minimum 
prompting  mode  requests  data  elements  by  displaying  the  data  element  mneumonics; 
while  the  maximum  prompting  mode  defines  each  of  the  data  elements. 

In  either  the  interactive  or  the  data  file  mode  the  IDP  creates  a 
formatted  input  file  for  the  SAASIII  analysis  program.  This  file  which  is 
stored  on  a mass  storage  device  must  be  identified  by  means  of  a file  name.  The 
program  thus  requests  the  user  to  enter  a name  for  the  file.  The  IDP  will  then 
catalog  the  file;  and,  an  error  will  result  if  the  file  already  exists.  The 
user  should  note  that  the  file  is  cataloged  as  a private  file  and  may  be  named 
using  read  keys  as  follows: 

QUAL I F I ER*NAME/ KEY1 / KEY2/ 

wheri  KEY1  is  the  write  key;  KEY2  is  the  read  key.  The  user  is  cautioned  to 
retain  these  keys  in  a secure  place,  if  used. 


7-2 


1 


7.2.2  Data  File  Mode 

In  this  mode  the  user  prepares  a data  file  containing  the  responses 
to  the  IDP  either  off-line  using  a paper  tape  or  tape  cassette;  or  on-line 
using  the  text  editor.  The  data  elements  are  described  by  Table  7.1.  The  pur- 
pose of  this  input  mode  Is  to  allow  the  user  to  prepare  an  Input  file  for  the 
IDP  so  that  the  prompting  produced  by  the  IDP  can  be  bypassed. 

The  IDP  will  query  the  user  to  determine  if  the  data  file  mode  is  to  be 
used.  If  so,  the  user  is  requested  to  enter  the  name  of  the  file  containing  the 
prepared  file  of  unformatted  input.  The  IDP  will  then  read  from  this  file  to 
satisfy  read  requests  completely  by-passing  the  Interactive  prompting. 

Occasionally  the  user  may  make  an  error  in  preparing  the  data  file 
which  will  cause  a read-error  in  the  IDP,  and  an  error  exit.  The  information 
printed  by  the  system  at  this  time  will  be  of  little  assistance  in  identifying 
the  data  elements  which  are  in  error.  A recommended  procedure  at  this  point  is 
to  place  the  IDP  in  the  prompting  mode;  but  to  cause  the  IDP  to  read  the  pre- 
created data  file.  This  can  be  accomplished  by  responding  "NO"  to  the  question 
"DO  YOU  WANT  TO  ENTER  DATA  FROM  A DATA  FILE?".  Then  in  response  to  the  first 
prompt  for  data  the  data  file  is  added  by  the  command 
mD  QUALIFIER*NAME 

The  IDP  now  prints  all  prompting,  but  reads  the  data  from  the  data  file. 
It  is  then  possible  to  identify  the  data  elements  in  error  by  noting  at  which 
point  the  error  exists. 

7.3  Examples 

The  SAASIII  examples  are  taken  from  reference  (11). 
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7.3.1 


Hollow  Cylinder  with  Uniform  Internal  and  External  Pressure 
(Lame1  Cylinder) 

The  four  element  idealization  of  the  cylinder  is  shown  by  Figure 


z 


Figure  7.1  Four  Element  Idealization  of  Cylinder 

The  loads  and  physical  parameters  are  as  given  below 
o.j  = 5000  psi 

oQ  = 10,000  psi 
E = 3 x 106  psi 
r^  ■ 1 in 

r ■ 2 in 
o 

v ■ .3 

7. 3. 1.1  Interactive  Mode 

The  interactive  Input  mode  for  the  creation  of  the  formatted 
SAASIII  input  file  is  shown  below,  where  the  formatted  file  is  to  be  called 
LAME 
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UOM  lies  3UIS9BI0 

ACCOUNT  NUMBER?  •BUI  SHAFER*  301  I4BHARRY.  SCHAEFFER.  MO 
PASSWORD? 


TAPE  USERS  SEE  *NEVS.F  OVERDUE-  TAPE* 

RUN  ID*  SHAPES  OS/M/74  10112133 

-tflIML  WMPOSIIM 
NEV  WORKSPACE  CREATED 
READY 

fXOT  ♦PREPROCESSOR. SAAS 

•44  INPUT  DATA  PREPROCESSOR  FOR  SAAS3  4*4 


INTER  NAME  OP  OUTPUT  DATA  PILE 
LAME 


DO  YOU  WANT  TO  ENTER  DATA  PORN  A DATA  PILE  CYES  OR  NO) 

JUL 

DO  YOU  WANT  MINIMUM  PROMPTING  (YES  OR  NO 


HOV  MANY  CASE ST 


WHAT  IS  THE  TITLE  FOR  THIS  CASE? 
LAME  CYLINDER 

•CONTROL  PARAMETERS- 


Ot TER  PLANE  STRESS/STRAIN  OPTION 

• BQ.O  NEITHER 
•EQ.t  PLANE  STRAIN 

• EQ.8  PLANE  STRESS 


ENTER  NUMBER  OP  NON-LINEAR  APPROXIMATIONS 
JL_ 

INTER  THE  MESH  GENERATION  PARAMETER 

• EG*  I FOR  MEM  GENERATION 

ENTER  THE  DATA  SET  < ITCARDS.  I PC  ARDS) 

• EQ.t  DATA  ON  NODE  CARDS 
•EQ.-l  DATA  ON  TAPE14 

• EQ.-e  CO) STAN T PROPERTIES 

INTER  DATA.. SET  ( IMAT.  I PRESS*  I SHEAR) 

IMAT  ■NUMBER  OP  DIFFERENT  MATERIALS 
I PRESS-NUMBER  OP  BOUNDARY  PRESSURE  CARDS 
1 Si  EAR-NUMBER  OP  BOUNDARY  MEAR  CARDS 

1.4.4 
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BITER  DATA  SET  ( TREF#  MTC#  NFREQ  > 

TWT  « REFERENCE  TEMPERATURE 

MTC  ■WMER  OP  TENSI ON /COMPRESSION  APPROX 

NFREQ  ( CALCULATE  FREQUENCY  IP  *EQ.t> 

IO.T 


•MESH  GENERATION  PARAMETERS- 
ENTER  THE  DATA  SET  I 

CMAXI#  MAX  J#  NLSEG*  NBC#  NMAT#  NIT#  PCI#  PCJ#  I CM.  JCM> 

MAXI  - MAX  I 
MAXJ  - MAX  J 

NLSEG  - NO  OP  LINE  SEGMENT  CARDS 
NBC  - NO  OP  BOUNDARY  CONDI  TITOS  CARDS 
NNAT  - NO  OP  MAT  BLOCK  CARDS 
NIT  - NO  OP  ITERATIONS 
PCI  • POLAR  PARAMETER  I 
PCJ  • POLAR  PARAMETER  J 
I CM  - I CURV  MOD 
JCN  - J CURV  MOD 
3#  3#  A#  B»l»0#0#0#0#0 

LINE  SEGMENT  1NPORMATION- 

POINT  COORDINATES  IN  THE  <1#J>-SPACE  AND  CR#Z>-SPACE  ARE  TO  BE 
ENTERED  ACCORDING  TO  THE  VALUE  OP  A CONTROL  PARAMETER#  N • 

IF  .EG.  0 ENTER  ONLY  FIRST  POINT  TO  DEFINE  A POINT 

IP  .EQ.  t ENTER  FIRST  AND  SECOND  TO  DEFINE  A LINE 

IP  #EQ#  8 ENTER  FIRST  AND  SECOND  TO  DEFINE  A DIAGONAL 

IP  .EG.  3 ENTER  THREE  POINTS  TO  DEFINE  CIRCULAR  ARC  WITH 
FIRST  AND  LAST  DEFINING  ENDS 

IP  .EG.  4 ENTER  <R#Z>  ONLY  TO  DEFINE  ARC  WITH  THIRD  PGINT  AS  CENTE 
R 

BITER  4 DATA  SETS 
CN> 

<II#JI#RI#Z1#...#I3#  J3#  R3#  Z 3> 
l#I.I..I.#l#3.  l.#2. 

I#  3# I »#8*» 3# 3# B«#  ft. 


3#l#t.#i.#l# !#!•#!• 


INFORMATION- 


DATA  IS  READ  TO  ASSIGN  BOUNDARY  CONDITION  TO  A BLOCK  OF  ELEMENTS 
BOUNDED  BY  II#I8#J1#J8. 

BITER  ..  8 DATA  SETS  « 1 1#  18#  Jt#  J8# CODE# XR# XE > 

I INI*  OR  JI«J8  TO  DEFINE  A LINE 
I|a!8  AND  J|aJ8  TO  DEFINE  A POINT 

XR  AND  XE  ARE  FORCE  AND/OR  DISPLACEMENTS  ACCORDING  TO  CODE  VALUE 
AS  FOLLOWS! 


CODE 

I 

t 

3 


XR 

R-LOAD 

R-DISP 

R-LOAD 

R-DISP 


XE 

E-LOAD 

Z-LOAD 

E-DISP 

E-MSP 


1*3#  1/  1*  8«  0«  * t» 

3«8^0»»0« 

•MATERIAL  BLOCK  INFOMATION- 

OITSR  I DATA  SETS  (MATI D#  MINI  #MAXI#MIMJ#NAXJ#  ALPHA) 

NAT 10  -MATERIAL  I.D. 

MINIM  MINIMUM  I 
MAX In  MAXIMIM  I 
MINJM  MINIMUM  J 
MAX J-  MAXIMUM  J 

ALPHAS  MATERIAL  PROPERTY  INCLINATION  ANGLE 

-material  Property  in format i cn- 

THE  PROGRAM  ALLOWS  FOR  THE  GENERAL  CASE  OP  A TEMPERATURE  VARYING 
ORTHOTROPIC  MATERIAL  WITH  DIFFERENT  PROPERTIES  IN  TENSION  AND 
COMPRESSION*  FOR  THE  GENERAL  CASE  THE  MATERIAL  PROPERTY  INFORMATION 
IS-  TO  BE  ENTERED  USING  FREE  FORMAT  IN  EXACTLY  THE  ORDER  SPECIFIED  IN 
THE  SAAS- 1 1 1 USERS  MANUAL!  WITH  EACH  LOGICAL  CARD  TERMINATED  BY  A 
CARRIAGE  RETURN* 

FOR  THE  CASE  OF  AN  ISOTROPIC  MATERIAL  HAVING  THE  SAME  PORPERTIES 
IN  TENS  KM  AND  CCMPRESSION  YOU  MAY  ELECT  TO  INPUT  A MINIMUM  OF  DATA 
DO  YOU  WISH  TO  ELECT  THE  LATER  OPTION? 

YES, 

FOR  EACH  OF  1 MATERIALS  ENTER  THE  FOLLOWING  DATA  SET 
CMIP#NT.  HHO#  TEX#  FOR) 

MID-  MATERIAL  I*D*  NUMBER 

NT  ■ NUMBER  OF  TEMPERATURES  AT  WHICH  NATL  PROPERTIES  ARE  DEFINED 
RHO»  NATL  DENSITY  <•  IF  5TATIC  PROBLEM) 

TEX-  THERMAL  EXPANSION  PARAMETER  <!F  I#  FREE  TNEfMAL  STRAIN# 

IF  0#  COEFFICIENT  OF  THERMAL  EXPANSION) 

POR-  EFFECTIVE  POROCITY 
FOLLOWED  BY  THE  DATA  SET 

<<  T(  I )#  E(  I )«.NU<  I )•  ALWAC  I )#  SIGNAC  I )•  RATIO<  I)  )#  I -I. NT) 

ALPHA  • FREE  THERMAL  STRAIN  OR  COEF  OF  THERMAL  EXPANSION 
FOR  TEX  • 1 OR  6!  RESPECTIVELY 
SIGMA  ■ YIELD  STRESS 

RATIO  - RATIO  OF  EFFECTIVE  PLASTIC  TO  ELASTIC  MODULUS 

illl  t/riltllQt  - 
j*#  3*  E6**3*fl»*Q**0* 

-BOUNDARY  PRESSURE  INFORMATION- 

ENTER  ..  4 DATA  SETS  C!#J#P> 

1 -NODE  I 
J -NODE  J 

P -NORMAL  PRESSURE 
7#4# SGOG* 

4# I#  S400* 

3#6#SGG0* 

6*9,10000* 

DO  YOU  WANT  A SUMMARY  OF  INPUT  DATA? 

NT 

XJBL 
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SUMMARY  OF  INPUT  DATA- 


NUMBER  OF  CASES  ■ I 


LAME  CYLINDER 

SOLUTION  OPTION  0 

PLANE  STRAIN  IP  1 
PLANE  STRESS  IF  8 

START  OPTIOI  I 

START  AT  BEGINNING  IF  1 
START  WITH  CON T CUR  PLOT  IF  8 
STOP  GPTICM  0 

AFTER  MESH  PLOT  IF  1 
BEFORE  CONTOUR  PLOT  IF  8 
DEFORMED  GRID  PLOT  IP  I 0 

PLOT  PARAMETER  0 

NO  OF  NOIL  INEAR  APPRO X 0 

MESH  GSMERATIOH  IF  1 1 

NO  OF  TEMPERATURE  CARDS  0 

NO  OF  NODAL  POINTS  0 

NO  CF  KLIMENTS  0 

NO  OP  PRESSURE  CARDS  0 

NO  OF  MATERIALS  I 

NO  OF  BOUNDARY  PRESSURE  CARDS  4 

NO  OP  BOUNDARY  SHEAR  CARDS  0 

REFERENCE  TEMPERATURE  Q. 

NO  OF  TB1./CCMP.  APPROX  G 

NATURAL  FREQ  PARAMETER  0 

-MESH  GENERATION  CONTROL- 


MAXI  a 3 

MAXJ  a 3 

NO  OF  LINE  SEQS  • 4 

NO  B*C*  CARDS  ■ . 8 

NO  OF  NATL  BLOCKS  .•  I 

NO  OF. ITERATIONS  a 0 

PCPI  a 0. 

PCPJ  a . §• 

I-CURV  a 0 

>CURV  a 0 

-LINE  SHBMENT  DATA- 


COORDS 

OF 

1ST  POINT 

COORDS 

CF  8MD 

POINT 

COORDS  OF 

3RD  POINT  PARAM 

I J 

R 

Z 

I J 

R 

Z 

I J R 

Z 

1 

1 

1.000 

1.000 

1 

3 

1.000 

8.000 

0 

0 

• 000 

• 000 

1 

1 

3 

1 .000 

8.000 

3 

3 

8.000 

8.000 

0 

0 

• 000 

• 000 

1 

3 

3 

8.000 

8.000 

3 

1 

8.000 

1.000 

0 

0 

• 000 

.000 

1 

3 

1 

8.000 

1.000 

1 

1 

1.000 

1.000 

0 

0 

• 000 

• 000 

1 

•BOUNDARY  CONDITION  DATA- 

LI  UI  LJ  UJ  CODE  R XS 


13  11  8.  0*  0* 

1 3 3 3 8.  0*  0* 


I 

I 

! 


•MATERIAL  BLOCK  AS  SI  (MM  ENT 
NATIB  LI  UI  LJ  UJ 


ALMA 


113  13  0* 

•MATERIAL  PROPERTY  INFORMATION- 

1HERE  IS  CARE>  I DIFFERENT  MATERIAL  GROUPS 


SATA  FOR  MATERIAL  NUMBER  1 FOLLOWS 

MATERIAL  I.D.  NUMBER— ——————— — 1 

NO  OP  TEMP.  AT  W1CN  PROPERTIES  ARE  DEFINED  1 

NASS  DSMSITT— — — — ..........  .60990 

THERMAL  EXPANSION  PARAMETER— .90090 

EFFECTIVE  POROSITY— —————  .00090 

MATERIAL  .SYMMETRY  PARAMETER 9 

TEMPERATURE  • .00000000 


TENSILE 

COMPRESSIVE 

THERM /YIELD 

E-N 

•30000*07 

.30000*07 

ALPKA-M 

•00000 

E-M 

•00000 

•00000 

ALPHA- N 

•00000 

- 

B-TH 

•09000 

•00000 

ALPHA- T 

• 00000 

NU-NN 

• 30000 

•30060 

SIONA-M 

•00000 

MU-MT 

•00000 

.00000 

SIQNA-N 

•00000 

MU-MT 

•00000 

.00000 

SIGMA- T 

• 00000 

e-NN 

•60000 

•00000 

PEMR 

•00000 

-NOUMDART  PRESSURE- 


MODE-  I MODE- J PRESSURE 

7 4 . 50000*04 

4 I .50000*04 

3 4 .50900*94 

6 9 .10000*05 

END  OP  CASE 
END  OP  DATA 


NORMAL  EXIT.  EXECUTION  TIME*  eOtS  MILLISECONDS. 


We  note,  turn  the  turn maiy  the  input  data  that  boundary 
p> itauAe  cL&io dated  with  point  (3,6)  i&  intellect;  the  value 
5 000.  should  be  10000.  We  w ill  theiefioie  edit  the  iile  LAME., 
and  change  the  value  o & pieAAuie. 


— AEP  uy». 

ED  Mb  09/00*  lfil  37-<00 J I F 
EDIT 

♦2  M 
l 

LANE  CYLINDER 


0 1 

• 0 0 

1 0 0 0 0 1 

4 

0 • 

00  0 

0 

3 

3 4 0 

1 0 *00000  *00000 

0 

0 

1 I 

1*M  1*00 

1 3 1*00  8.00  0 

0 

• 000 

• 000 

1 

1 3 

i*oo  e*oo 

3 3 8*00  8*00  0 

0 

• 000 

• 000 

1 

3 3 

8*00  0*00 

3 1 8*00  1.00  0 

0 

• 000 

• 000 

1 

3 1 

0*00  1*00 

1 t 1*00  1.00  0 

0 

• 000 

• 000 

1 

1 

3 1 ! 

8*0000  *00000  .00000 

1 

3 3 3 

8*0000  0 00800  *00000 

1 

1 3 1 

3 *000 

1 

1 *00000 

•00000  *00000  8 

•00000  • 30000407 

•00000  *80000  *30000 

•00000 

•00000 

• 

•0000 

* 

• 30000407 

•00000  *00000  *30000 

•00000 

•00000 

• 

•0000 

• 00000 

•00000  *00000  *00000 

•00000 

•00000 

• 

•MOO 

7 4 5000*0 

4 I 5000*0 

3 0 5000*0 

0 9 10000* 

END  OF  CASE 
END  OF  DATA 
EOT  AT  LINE  01 
• 


We  can  now  (,-tnd  and  change  the  appropriate  tine  LAME. 
by  meam  o&  the  fiolCowing  edit  aommandA 


4 

3 6 5000*0 

»C  75/10/ 

3 6 10000*0 

• E 

■D  EDIT  01  LINES  OUTPUT 


7. 3. 1.2  Data  File  Mode 

The  input  data  for  the  SAASIII  will  be  created  in  the  element 
LAME  of  the  file  SAASPROBS  by  using  the  text  editor.  First  a file  SAASPROBS 
must  be  created  by  using  the  system  commands: 

SCAT,  UP  SAASPROBS. 

OFREE  SAASPROBS. 

OASG.A  SAASPROBS. 

The  data  is  now  input  by  means  of  the  text  editor  as  follows. 

<?ED  SAASPROBS.  LAME 
INPUT 

(The  input  data  is  now  shown) 

The  IDP  is  now  executed  in  the  data  file  mode  as  follows,  where, 
again  we  assume  that  LAME  contains  the  formatted  output  file. 


•XttT  ♦PREPROCESSOR. SAAS 


•••  INPUT  DATA  PREPROCESSOR  POR  SAAS3  ••• 


MTER  NAME  OP  OUTPUT  DATA  FILE 
LAME 


DO  YOU  WANT  TO  ENTER  DATA  FOFM  A DATA  FILE  (YES  OR  NO> 


■ITER 


ft  OP  INPUT  BATA  FILE 


1 


- SUMMARY  OP  INPUT  BATA 
NUMBER  OP  CASES  • 


LAME  CYLINDER 

SOLUTION  OPTION  0 

PLANE  STRAIN  IP  1 
PLANE  STRESS  IP  & 

START  OPTI ON  1 

START  AT  BEQINNIN6  IP  t 
START  WITH  CONTOUR  PLOT  IP  2 
STOP  OPTION  0 

AFTER  MESH  PLOT  IP  t 
BEFORE  CONTOUR  PLOT  IF  2 
DEFORMED  GRID  PLOT  IF  1 0 

PLOT  PARAMETER  0 

NO  OF  NONLINEAR  APPROX  • 

MESH  GENERATION  IP  1 1 

NO  OP  TMPERATURE  CARDS  0 

NO  QF  NODAL  POINTS  0 

NO  OF  ELEMENTS  • 

NO  OP  PRESSURE  CARDS  0 

NO  OP  MATERIALS  1 

NO  OP  BOUNDARY  PRESSURE  CARDS  4 

NO  OF  BOUNDARY  SHEAR  CARDS  0 

REFERENCE  TEMPERATURE  •• 

NO  OF  TIM./C(MP.  APPRO K 0 

NATURAL  FREQ  PARAMETER  0 

•MESH  GENERATION  CCMTROL- 


MAXI  ■ 3 

MAXJ  - 3 

NO  OF  LINE  SEGS  • 4 

NO  B.C.  CARDS'  • . 2 

NO  OF  NAIL  BLOCKS  .*  I 

MO  OF. ITERATIONS  - 0 

PC PI  * •. 

PCPU  • ..  O* 

I-curv  n t 

J-CURV*  0 

•LINE  S BEN  BIT  DATA* 


COORDS 

OF  1ST 

POINT 

COORDS 

OF  2ND  POINT 

COORDS 

OF  3RD 

POINT  1 

PARAM 

I J 

R 

C 

I J 

R 

z 

I 

J 

R 

Z 

1 1 

I.OfO 

1.000 

1 3 

1.000 

8.  000 

0 

0 

• 000 

• 000 

1 

1 3 

1.000 

8.000 

3 3 

8.000 

t.  000 

0 

0 

• 000 

• 000 

1 

3 3 

2.000 

8.000 

3 1 

2.000 

1.000 

0 

0 

• 000 

• 000 

1 

3 1 

8.000 

1 .000 

1 1 

1.000 

1.000 

0 

0 

• 000 

• 000 

1 

•BOUNDARY  CONDITION  DATA- 

LI  Ul  LJ 

UJ 

CODE 

R 

xz 

1 

3 1 

1 

2. 

0. 

0. 

1 

3 3 

3 

8. 

0. 

o. 
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-MATERIAL  UCK  ASSI6NMI 
■ATI  0 LI  01  LJ  UJ 


T< 


M.FHA 


113  13  0. 

-MAT  DUAL  PROPERTY  IN  FORM  ATI  CM- 

THERE  IS  (ARE)  1 DIFFERENT  MATERIAL  GROUPS 


DATA  FOR  MATERIAL  HUMBER  1 FOLLOWS 

MATERIAL  I.D.  NUMBER- 

NO  OP  TEMP.  AT  WICK  PROPERTIES  ARE  DEFINED 


MASS  DENSITY 


THERMAL  EX  PAN  SI  CM  PARAMETER- 

EFFECTIVE  POROSITY — 

MATERIAL  .SYMMETRY  PARAMETER- 


TEMPERATURE  • .00000000 


I 

1 

•setos 

•OOOBO 

•00000 

e 


TENSILE  OPPRESSIVE  THERM/YIELD 


E-N 

• 30000*07 

• 30000*07 

ALPHA  41 

•00000 

t 

E-M 

•00000 

•00000 

ALPKA-N 

•00000 

1 

E-TH 

•00000 

•00000 

ALPHA- T 

•00000 

» 

NU-MN 

•30000  t 

• 30000 

SIGMA-M 

•00000 

NU-MT 

•00000 

•00000 

SIGMA-N 

•00000 

I 

NU-NT 

•00000 

•00000 

SIGMA-  T 

•ooooo 

1 

t 

G-MN 

•00000 

•ooooo 

PEMR 

•ooooo 

-BOUNDARY  PRESSURB- 

I 

NODE- I NODE-J  PRESSURE 

i 7 4 • 50000*04 

4 1 • 50000*04 

3 6 • 10000*05 

6 9 .1 0000* OS 

BID  OF  CASE 
END  OF  DATA 


NORMAL  EXIT.  EXECUTION  TIME*  0149  MILLISECONDS. 

4 


I 


| 


I 
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7. 3. 1.3  Analysis  Program  Execution 


Using  one  of  the  two  modes  in  the  IDP  we  have  created  a file 
LAME,  which  contains  the  formatted  input  data  for  SAASIII.  The  analysis 
program  requires  more  core  than  is  normally  available  for  demand  execution. 
The  program  will  thus  be  sent  to  the  batch  queue  for  execution  by  means  of 
the  OSTART  command.  First,  however,  we  must  create  a run  stream  to  execute 
the  program.  Let's  call  this  run  stream  START-SAAS  we  then  proceed  as 
follows: 

(SCAT, UP  START-SAAS. 

@FREE  START-SAAS. 

@ASG,A  START-SAAS. 

The  output  of  the  SAASIII  program  will  be  copied  to  SAASOUT  using 
the  following  run  stream 

@ED  START-SAAS. 

INPUT 

#RUN  runid  , acct  no  , SYMPOSIUM,  5,  500 
#CAT,U  SAASOUT 
#FREE  SAASOUT. 

#ASG,A  SAASOUT. 

IBRKPT  PRINT3/SAAS0UT 
#XQT  PROGRAM. SAAS 
#ADD  LAME. 

#FIN 

The  #'s  are  now  changed  to  @'s  by  using  the  change  feature  of  the 

editor. 
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After  satisfying  yourself  that  the  file  is  correct  it  is  submitted 


I 

I 

T 


to  the  batch  queue  by  the  conmand 
OSTART  START-SAAS. 

7. 3. 1.4  Analysis  program  Output  Listing 

The  listing  of  the  results  of  SAASIII  for  the  Lame'  cylinder 
problem  is  as  follows. 


L : 


I 

i 


r 


CF  CE  KE'ATICN  IfT  C r J-  A"I  CN 


**  ajhvun 

VALUE 

r*  c-  - ^ * 

"nr  vrrK 

3 

MA  X * M’J'*! 

VAL’J" 

or  j in 

IMF  "EEH 

3 

M y m rrn 

CF  LINE 

etke  n 

7 CA~EE- 

4 

\»y  •.•'*  r -> 

or  "dut: 

0 A ~ Y CO 

: % - - I «.  Jll  l-n.K 

2 

•Ctf'E"  CF  VA"ZOIAL  SLCCK  CAROS- 


1 

0 


"OLA"  "CC~0  TN A"E  F AF.AN.ETE  F I - 

"OLA?  CC00''T\'AT2  " ARAMETE?  J 

” CURVATURE  VC^ir-CATZCN f 

J CUR VATU RE  NCEI'ICA'ICN C 


.ccee 
. coon 


- ‘try  ■»  t j 

vl  01 

"1 

12 

,r "»  <7  n 

c ~ 

7 r»  T 7 

4-  <-  A _ 

J7 

57  7 -* 

1 

i l.rco 

1.2 

CO  1 

1.2C" 

0 p -»n 

..  • L>  u v. 

o 

u 

« 1.  • 

CT=  ".  rL 

= r.  Errr 

= 2. 

rinc= 

.EEC  zr 

NO 

— f *■>  r> 

— ' • w U 

ITER: 

1 II NCr 

- i 

A? 

« -* 

1 *’ 

« 

i.cee 

IMF  U ' T1 

vl  FI 

’i 

12 

J2  R F 

Z 2 13 

J3 

• T 7 

• » . Z. 

l 

7 i.cco 

2.0 

CC  3 

7 2. EC" 

-»  . Q p r>  r, 

0 

• r c ° • ? 

C"-  2.  \ 

: c . err  ~ 

— o 

— ^ • 

RIf,Cr 

.EEC  zr 

NC 

~ .EEC 

ITER: 

l'  ~ X Z N C — 

: j 

A? 

AZ 

r*  t 

1 .7?" 

2.0CC 

INPUT  71 

v’l  FI 

21 

12 

J2  R 2 

Z2  13 

J3 

R 2 7 7 

3 

7 2. COO 

2.000  3 

1 2 .00° 

1.000  o 

0 

iOOO  .0 

OT-  0.  CJ 

= 2.  C IF  r 

=■  2. 

RINCr 

.CCO  ZI 

NC 

: - ■ 5C0 

ITER: 

1 I INC: 

I J 

A? 

AZ 

3 7 

2.nGr 

i.fec 

I*FUT  71 

J1  R1 

21 

12 

J2  R? 

Z2  I? 

J3 

R2  7 7 

3 

1 2. OCT 

1.? 

00  1 

1 1.000 

1.000  n 

0 

.000  .0 

DT  = 2.  Ev’ 

r r.  DIFF 

= 2. 

RIRC  = 

-.EEC  ZI 

NC 

: .ore 

ITER: 

1 I INC:  - 

7 J 

A" 

AZ 

2 1 

l.rCC 

1.0CC 

* 

, 

c^creinatee 

CALCULATE!? 

AFTER 

43  ITER 

ATICNS 
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TABLE  7.1  SAASIII  Preprocessor  Input  Data  Elements 


1.  Title  (.LE.72  Characters) 

2.  Plane  stress/straln  option 

3.  Number  of  non-linear  approximations 

4.  Mesh  generation  parameter  (MGP) 

5.  (NODES,  NELMTS) 

6.  (ITCARDS,  IPCARDS) 

7.  (IMAT,  IPRESS,  ISHEAR) 

8.  (TREF,  NTC,  NFREQ) 

9.  IF  (NFREQ. EQ.O)  enter  no  data,  otherwise  enter  (R,Z) 

10.  IF  (MGP.NE.l)  enter  no  data,  otherwise  enter  the  data  set 

10. A (MAXI,  MAXJ,  NLSEG,NBC,NMAT,NIT,PCI ,PCJ .ICM.JCM) 

Enter  NLSEG  sets  as  follows 

10.B.1  N 

10. B. 2 (I|i  , R-j , Zj , ... , 1^,  Jj,  Rj,  Zj) 

Enter  NBC  data  sets  as  follows 
10. C (II,  12,  Jl,  02,  ICODE,  XR,  XZ) 

Enter  NMAT  data  sets  as  follows 

10.  D (MATID,  MINI,  MAXI,  MINJ,  MAXJ,  ALPHA) 

11.  IF(ITCARD.EQ.O)  enter  no  data  otherwise 

ll.A.l  if  (ITCARD.GT.O)  enter  ITCARD  sets  (R,Z,T) 

11. A.2  If  ( ITCARD. EQ. -2)  enter  (T) 

12.  IF  (IPCARD.EQ.O)  enter  no  data,  otherwise 

12. A. 1 if  (IPCARD.GT.O)  enter  IPCARD  sets)  (R,Z,P) 

12. A. 2 If  (IPCARD. EQ. -2)  enter  P 

IF  (MGP.NE.l)  proceed  to  Item  15 

7-26 
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TABLE  7.1  (Continued) 


13.  Enter  NODES  date  sets  (NODE,  CODE,  R,  Z,  XR,  XZ,  T (If  ITCARD.EQ.Q) , 
PC  If  IPCARD.EQ.O)). 

14.  Answer  query  "Do  you  Input  minimum  data?" 

14. A YES  Enter  IMAT  data  sets 

(MID,  NT,  RHO,  TEX,  POR) 

(Tr  Ej , Nil j , ALPHA1 , SIGMA^  RATIO^  1 » 1 , NT 

14.  A NO 

Enter  IMAT  data  sets 

(MID,  NT,  RHO,  TEX,  POR) 

followed  by  NT  data  sets  for  each  MID: 

(T’  EMt*  Ent’  Eet*  vMNt’  wMet’  vNet‘  E45t) 

( eMc’  ENc’  Eec*  vMNc*  vMet*  vNec’  E45cJ 
( oM,  oN,  a0,  aM,  oN,  aQ) 

15.  enter  (IPRESS)  data  sets 

(I,  J,  PRESS) 

16.  enter  (ISHEAR)  data  sets 

(I,  J,  SHEAR) 

17.  Answer  to  query  "Do  you  want  sunroary  of  input  data?" 
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